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ABSTRACT 

Context. There is gathering indirect evidence suggesting non-conservative evolutions in Algols. However, the systemic mass-loss rate 
is poorly constrained by observations and generally set as a free parameter in binary-star evolution simulations. Moreover, systemic 
mass loss may lead to observational signatures that are still to be found. 

Aims. Within the ‘hotspot’ ejection mechanism, some of the material that is initially transferred from the companion star via an 
accretion stream is expelled from the system due to the radiative energy released on the gainer’s surface by the impacting material. 
The objective of this paper is to retrieve observable quantities from this process, and to compare them with observations. 

Methods. We investigate the impact of the outflowing gas and the possible presence of dust grains on the spectral energy distribution 
(SED). We used the ID plasma code Cloudy and compared the results with the 3D Monte-Carlo radiative transfer code Skirt for 
dusty simulations. The circumbinary mass-distribution and binary parameters are computed with state-of-the-art binary calculations 
done with the Binstar evolution code. 

Results. The outflowing material reduces the continuum flux-level of the stellar SED in the optical and UV. Due to the time- 
dependence of this effect, it may help to distinguish between different ejection mechanisms. Dust, if present, leads to observable 
infrared excesses even with low dust-to-gas ratios and traces the cold material at large distances from the star. By searching for such 
dust emission in the WISE catalogue, we found a small number of Algols showing infrared excesses, among which the two rather 
surprising objects SX Aur and CZ Vel. We find that some binary B[e] stars show the same strong Balmer continuum as we predict 
with our models. However, direct evidence of systemic mass loss is probably not observable in genuine Algols, since these systems 
no longer eject mass through the hotspot mechanism. Furthermore, owing to its high velocity, the outflowing material dissipates in a 
few hundred years. If hot enough, the hotspot may produce highly ionised species such as Siiv and observable characteristics that are 
typical of W Ser systems. 

Conclusions. Systemic mass loss, if present, leads to clear observational imprints. These signatures are not to be found in genuine 
Algols but in the closely related p Lyrass, W Serpentis stars. Double Periodic Variables, Symbiotic Algols and binary B[e] stars. We 
emphasise the need for further observations of such objects where systemic mass loss is most likely to occur. 

Key words. Binaries: general - Radiative transfer - Stars: circumstellar matter - Methods: numerical 


1. Introduction 

Algol systems are short-period, semi-detached, binary systems 
(with orbital periods from several hours to tens of days) gen¬ 
erally composed of a hot (B-A) main-sequence star and a less 
massive but more evolved companion overfilling its Roche lobe 
and transferring mass via Roche lobe overflow (RLOF; Kopal 
1955;Giuricinetal. 1983). 

Although Algol systems have been studied since the 18* 
century (Goodricke 1783) and seem to be rather simple objects 
within the zoo of binary systems, it is still unclear whether the 
mass transfer can be fully conservative, or if systemic mass loss 
is a universal feature of Algols as it has been inferred already by, 
e. g., Crawford (1955). It has been shown that systemic angular- 
momentum losses (and thus mass losses) must be involved to 
reproduce the current state of some observed Algols. For ex¬ 
ample, Refsdal et al. (1974) found that only a non-conservative 
evolution can reproduce the system AS Eridani. Massevitch 
& Yungelson (1975) noticed that Algol models must loose at 
least 40% of their mass in order to reproduce observed proper¬ 
ties. Eventually, Sarna (1993) concluded that the Algol proto- 
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type p Per must have lost 15% of its initial total mass during its 
evolution and 30% of its total angular momentum. Statistically, 
Chaubey (1979) found that semi-detached binaries have lower 
orbital angular momentum than detached systems of the same 
total mass. Finally, many other clues (e. g., comparison of grid 
of models with observations or population synthesis analysis; 
Mezzetti et al. 1980; Nelson & Eggleton 2001; van Rensbergen 
et al. 2011 ) converge in favour of non-conservative evolution. 

Different mechanisms have been proposed to explain sys¬ 
temic mass loss. They may be classified into four groups: bipolar 
jets (Umana et al. 2000, 2002; Ak et al. 2007; Lomax & Hoff¬ 
man 2011), enhanced winds (chapter 4 of Maeder 2009; Tout 
& Pringle 1992), losses through the outer Lagrangian point -C 3 
(Sytov et al. 2007), and ‘hotspot’. The latter mechanism appears 
to be the most promising, because the mass-loss rates derived 
from the first three mechanisms are too low to account for the 
inferred systemic mass loss of some Algols. In this scenario, the 
mass loss is driven by the radiation pressure of a hotspot (van 
Rensbergen et al. 2008, 2010; Deschamps et al. 2013, hereafter 
Paper I) that forms when the stream of material coming from the 
donor via RLOF hits the accreting star and releases energy to 
heat the stellar surface and increase the luminosity locally. This 


Article number, page 1 of 21 




0.25 


Roche lobe 



Fig. 1. Schematic overview of an Algol system undergoing systemic 
mass loss due to a hotspot on the surface of the gainer. The ticks and 
numbers indicate the different phases (0.0: primary eclipse). 


hotspot is observed in some Algol-like stars (e. g., Banks et al. 
1990; Welland et al. 1995; Richards & Albright 1999, in case of 
direct accretion on the stellar surface and Harmanec 2002; Soy- 
dugan et al. 2007 in case of accretion on the edge of a disc). 

In the line of van Rensbergen et al. (2008), Deschamps et al. 
(2013) re-derived a model for systemic mass loss via hotspot 
based on energy conservation at the impact location and showed 
that for a typical initial 3.6 Mq + 6 Mq system with an ini¬ 
tial orbital period P = 2^5, a large amount of matter (of the 
order of 1.4 M©, i. e. 15 % of the initial total mass of the sys¬ 
tem) can be lost during the short-lasting (2 x 10^ years) rapid 
(44rlof = 1X10“^-1 X10“"^ M© yr“') mass-transfer phase. This 
total amount of mass lost is in good agreement with the value 
inferred by Sarna (1993) to reconcile Algol models with obser¬ 
vations. Using the mass-loss history, and assuming the wind 
to follow the prescription of Castor et al. (1975) for line-driven 
winds, a density profile for the material surrounding the system 
was derived. 

In this paper, we confront this prediction of systemic mass 
loss to observations by converting the density profile into ob¬ 
servable quantities. We will thus perform an analysis of light 
emission and absorption by this material using the codes Cloudy 
(Ferland et al. 2013) and Skirt (Baes et al. 2003,2011 ; Camps & 
Baes 2015). Section 2 briefly recaps and discusses the assump¬ 
tions made to obtain the density profile and presents the char¬ 
acteristics of the model system and the set-up of the two codes. 
The results are presented in Sect. 3. We then confront our find¬ 
ings to Algol observations and broaden our study to promising 
Algol-related objects in Sect. 4 before concluding in Sect. 5. 


2. The model system 

2.1. Density profile of the out-flowing matter 

We follow the ballistic trajectory (RLOF-stream) and thus the 
specific kinetic energy of the material escaping from the donor 
star via Roche-lobe overflow through the Lagrangian point Xi 
and accreting onto the gainer star by direct impact (Paper I). 
When impacting onto the stellar surface, a fraction of this en¬ 
ergy heats up the perturbed stellar surface, forming a hotspot 
(see Fig. 1), and the remainder is radiated away. Mass loss is 
assumed to occur whenever the total luminosity at the impact lo¬ 
cation (stellar -H hotspot luminosity) exceeds the local Eddington 


luminosity. The latter is computed using the gainer upper layers’ 
opacity. 

We assume the wind (hereafter outflow) emitted by the 
hotspot to be line-driven, and follow the Castor et al. (1975) 
prescription to compute the outflow speed iiflow Due to the or¬ 
bital motion, the outflow will exhibit a spiral shape (Fig. 1; see 
Sect. 2.2). The acceleration of the particles moving in this spi¬ 
ral outflow will be considerably reduced once the matter leaves 
the “irradiation” cone defined by the opening angle 6(flow (whose 
estimate is provided in Appendix A). It is important to stress 
that a mass element ejected from the hotspot will escape the sys¬ 
tem, even if the outflowing material is not constantly lit up by 
the hotspot, because it exceeds its escape velocity after a small 
fraction of the orbital period' (0.038 x Porb at the beginning 
of the mass-transfer phase, 0.011 x Porb at the end of the non¬ 
conservative phase), during which the change of the orientation 
of the hotspot is negligible. 

The density of the material surrounding the system at a given 
time t and distance r is calculated using the continuity equation 
(see Paper I for a complete description): 


p(r, t) = 


4Tflow(^i 0 
Anf^r^vno^ir, t) ’ 


( 1 ) 


where Mflow is the systemic mass-loss rate calculated with the 
Binstar code (Siess et al. 2013, Paper I) and 4-nfc is the solid 
angle into which the material is ejected, /c is the covering factor 
defined by Eq. A.9 and is related to the opening angle Ohow of the 
outflow (see Eig. 1). The possibility that the outflowing matter 
be blocked by the infalling material from the RLOE stream has 
not been considered, because the hotspot will be elongated as 
a consequence of stellar rotation (the accreting star is spun up 
and no longer synchronised, as observed by Mukherjee et al. 
1996). It is therefore very likely that there is only a small overlap 
between the stream and the outflow. 

By computing the density profile of the circumstellar ma¬ 
terial at different stages during the evolution, we can estimate 
the total amount of mass potentially detectable. Eigure 2 shows 
the evolution of this quantity during the non-conservative phase. 
We see that the total amount of detectable mass surrounding the 
binary is maximal before the reversal of the mass ratio (mass ra¬ 
tio ^ = 1 at f ~ 75 000 yr: the beginning of the Algol phase). 
The gainer surface temperature mainly evolves during the rapid 
phase of mass transfer and reaches a plateau when the mass ratio 
is reversed. 

Eor our radiative-transfer simulations, we selected three con¬ 
figurations (see Table 1): (i) right after the mass-ratio reversal 
(f = 80 000 yr; hereafter model A); (ii) at f = 120 000 yr (model 
B); (iii) at the end of the non-conservative phase {t - 162 000 yr; 
model C, our default configuration). In the density profiles of 
models A and B, the densities in the ionised region close to the 
system become so high that an increasing amount of the starlight 
is transformed to recombination lines and continua that are not 
explicitly calculated by Cloudy. Indeed, the model atoms used 
by Cloudy account for a finite number of energy levels and emis¬ 
sion from recombinations to higher levels cannot be included in 
the spectrum. As a result, a significant fraction of the energy 
released in these recombinations and in the subsequent cascades 
to the corresponding ground states are missing in the emergent 
spectrum. In these conditions, it is difficult to predict observable 
quantities, and we decided to use the density profile of the model 
C for the models A and B while keeping the hotspot and stellar 

’ As presented in Fig 3, the escape velocity is reached at a fixed dis¬ 
tance from the star and therefore after a fixed amount of time. 
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Fig. 2. Red crosses: evolution of the mass surrounding the system 
(integrated up to the radius where the density of the outflow equals the 
ISM density; see text). The cyan circles (labelled with the model name) 
indicate the models selected for the simulations presented in this work. 
Long-dashed green line (right hand scale): evolution of the gainer-star 
surface temperature. Magenta dotted line and dashed blue line: evolu¬ 
tion of the masses of the donor and gainer (in unit of Mo/500). 

parameters given by the evolutionary model. Therefore, the sim¬ 
ulations corresponding to cases A and B should be considered 
as exploring the sensitivity of our results to the luminosity of 
the hotspot, which is the property differing the most between 
cases A, B and C (Table 1). The sensitivity to the circumstellar 
shell mass (which was the initial goal for considering cases A, 
B and C - but which could not be fulfilled for the reason ex¬ 
plained above) may instead be probed by varying the dust-to-gas 
ratio O. As discussed in Sect. 2.3, a simulation for case C has 
been run with 200 times more dust than standard (model labelled 
’C_hd_0.7r in Table 2). This case thus mimicks a situation with 
a thicker dust shell. 

Figure 3 shows the density profile of the outflow for differ¬ 
ent covering factors (/c) for model C. The profile is cut at the 
distance r^ut, where the density drops below the ISM value. This 
occurs at 6.08 x lO'^ cm (» 40600 AU or 2.5x10^ Aorb, where 
Aoib is the orbital separation) from the gainer for /c = 0.71. We 
do not take into account any bow-shock formation at this loca¬ 
tion. Note that when the systemic mass loss stops, within 300 
years the material disperses and the density profile drops every¬ 
where below the ISM density 10“^^ g cm“^) because of the 
very high velocity of the line-driven outflow (Fig. 3). 


Fig. 3. Density (left-hand scale, solid lines) and velocity (right-hand 
scale, dashed line) profiles of the outflowing material (similar - but 
with different parameters - to Fig. 9 of Paper I). Blue line: spherically- 
symmetric and isotropic outflow corresponding to a covering factor of 
/c = 1 (see text. Sect. 2.1); green lines: ejection conhned in a disc of 
opening angle 90° (/c = 0.71); red line: ejection confined in a disc of 
opening angle 74° (/c = 0.6); The solid light-grey lines are the density 
profiles for = 0.71 at 100 different times uniformly distributed during 
the non-conservative phase. 


the hotspot and will decrease as the orbital motion moves the 
hotspot out of sight. Ultimately, right before the hotspot points 
again toward the observer, the density close to the star reaches a 
minimum as only the material that has been launched orbital pe¬ 
riods ago is seen far away from the star. To mimic this effect and 
investigate any variability due to the orbital motion, we perform 
some simulations with a reduced, constant density profile close 
to the star (in the following, we call these simulations ‘plateau’ 
and label them with a ‘p’ in the name). Considering the orbital 
motion (Porb = 8^33 at the end of the non-conservative phase), 
we can estimate the density of the outflow ahead of the hotspot. 
Given the travel time (one orbital period) and the flow veloc¬ 
ity, we can estimate the distance travelled by the material, rpi, 
as well as the density at that location, ppi. As a first approxi¬ 
mation, we adopt a density profile where the density is equal to 
Ppi up to Tpi and then follow the standard (unmodified) density 
profile. For the model with = 0.71, we retrieve a radius of 
Tpi = 6.05 X 10'^ cm (si 4.05 AU or 25 Aorb) and a corresponding 
density of ppi = 9.44 x 10^*^ gcm“^. 


2.2. Spiral structure 

Due to the orbital motion and the fact that the outflow comes 
from a small area on the star, we must expect a spiral structure 
to form close to the star. The width of the spiral arm increases as 
a function of the distance from the star (Fig. 1 ) and eventually the 
spiral arm will overlap itself so that the spiral structure fades to 
give rise to a disc instead. Using an opening angle for the outflow 
of 90° (standard 0flow configuration corresponding to /c = 0.71) 
and the velocity profile shown in Fig. 3, we estimate the distance 
at which the spiral fades to be » 3.6 AU (22 Aorb)- This dis¬ 
tance is reached by the outflow on the time scale of an orbital 
period. At larger distances, our assumption of a smooth density 
profile is thus sufficiently accurate. At smaller distances, we ex¬ 
pect deviations due to the spiral structure. In these regions, the 
density will be the highest when the line of sight points towards 


2.3. Chemistry of the outflow: gas composition and dust 
formation 

One important point for the radiative-transfer analysis is the 
composition of the gas, which we assume to be that of the donor 
as listed in Table D. 1 . The main signature is the depletion of car¬ 
bon and enhancement of nitrogen during the non-conservative 
phase (since the mass transfer reveals the inner layers of the gi¬ 
ant star). 

Although grains can be useful tracers of the cool outermost 
regions of the outflow, it is not clear whether or not they can 
form. Indeed, the escaping material cools as it leaves the gainer’s 
vicinity. Eventually, the temperature drops below the condensa¬ 
tion temperature, allowing for the formation of grains. However, 
UV and X-ray photons from the hotspot and star (e. g., Ignace 
et al. 2008) may prevent grain formation. Surprisingly, grains 
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are found in the surroundings of Wolf-Rayet stars (Barniske et al. 
2008), which are also emitters of hard radiation. 

As there are no observational constraints on the chemical 
composition and optical properties of (hypothetical) dust grains 
around Algol systems, we use the standard ISM model imple¬ 
mented in Cloudy. However, we exclude carbonaceous grains, 
because the evolutionary stages of both stars imply C/O < 1 (see 
Table D.l) so that we can assume that most carbon atoms in the 
outflow are locked up in CO molecules. For the slope of the size 
distribution of the remaining silicate grains we use the expres¬ 
sion suggested by Mathis et al. (1977), 


dUgr 

dfl 


oc a 


3.5 


( 2 ) 


where Ugr is the number density of grains and a the grain size. 
The size distribution is discretised in ten bins, which is the maxi¬ 
mum number of size bins allowed by Cloudy, ranging from 0.005 
to 0.25 pm, with a ratio amax/flmin ~ 1.48 in each bin, where flmax 

and flmin are the upper and lower size limit of the bins, respec¬ 
tively. Within each bin, the optical properties of the dust grains 
are averaged over the respective size interval. For our simula¬ 
tions, the chosen number of bins and their width and location 
have very little impact on the results. Varying the number and 
size of the bins only results in a deviation of ^ 2% in the flux ra¬ 
tios and magnitudes. Barniske et al. (2008) evaluate the impact 
of the lower cut-off on the grain size (0.005 pm in our case) and 
And that only the flux blueward of 8.3 pm is altered, which has 
no influence on WISE bands 3 and 4 (> 10 pm) that will be used 
in Sect. 4.2.2. Similarly, we ran simulations with different grain 
sizes and found that the region of the spectra below ~6 pm does 
not depend on grain size so that 2MASS fluxes (< 5 pm) are not 
affected. 

The optical properties of the grains are calculated by use of 
Mie theory as described by van Hoof et al. (2004). As the gainer 
star is of spectral type B, its radiation will probably interfere with 
the formation of dust grains, and we cannot expect the dust-to¬ 
gas ratio (O) to be as high as in the ISM. Instead, we use the 
results of Barniske et al. (2008), who And dust in the surround¬ 
ings of WN stars. From their best-fit models we derive values of 
<5 up to 2 X 10“^, which we adopt in this work for our standard 
dusty configuration (models labelled ‘Id’ for low-dust). To in¬ 
vestigate the impact of the dust-to-gas ratio, we also computed a 
model with O = 4 x 10“^ (model labelled ‘hd’ for high-dust). 

We do not explicitly calculate the destruction of dust grains 
caused by the high temperatures close to the star. Instead, we 
use the Cloudy code to derive a temperature profile for a pure- 
gas model (models labelled ‘df’ for dust-free) and then define a 
dust-free region around the star where the temperature exceeds 
the dust sublimation temperature (^ 1400 K for silicate dust, see 
Martin & Rouleau 1991; Laor & Draine 1993). The circumstel- 
lar material will be dust-free up to a transition radius (radius 
at which the temperature goes below the condensation tempera¬ 
ture). 


2.4. Stellar parameters 
2.4.1. Luminosities 

The system is composed of the donor star (a late G-type red gi¬ 
ant), the gainer star (a B-type dwarf), the hotspot and the material 
escaping from the hotspot (Fig. 1). The properties of the system 
(masses, radii, stellar and hotspot luminosities, stellar temper¬ 
atures; Deschamps et al. 2013) are computed with the Binstar 
code and given in Table 1 for the three selected models. 


Table 1. Parameters for the three models; Model A: mass ratio reversal; 
Model B: intermediate case; Model C: end of non-conservative mass 
transfer (for this case, three hotspot temperatures are investigated). 


Model A 

Gainer Donor 

Hotspot 

Radius (R©) 

3.87 5.73 

0.14(0.19) 

Mass (Mq) 

4.65 4.26 

— 

Luminosity (L©) 2 496 54 

7 

Temperature (K) 20 840 6 555 

35 000 

Porb 

2.3 d 


■^orb 

15.5 Ro 



Model B 

Gainer 

Donor 

Hotspot 

Radius (R©) 

3.13 

6.19 

0.24 (0.26) 

Mass (Mq) 

5.66 

2.83 

— 

Luminosity (L©) 

1924 

38 

20 

Temperature (K) 21 780 

5 784 

35 000 

Porb 


3.3 d 


■^orb 


19.0 Ro 



Model C 

Gainer 

Donor 

Hotspot 

Radius (Rq) 

3.22 

9.55 

0.87 (0.63) (C XX 0.71 cold) 
0.64 (0.63) (C XX 0.71 inter) 
0.49 (0.63) (C_xx_0.71_hot) 

Mass (Mq) 

6.50 

1.68 

_ 

Luminosity (Lq) 

2195 

63 

137 

Temperature (K) 

22150 

5 297 

30000 (C_xx_0.71_cold) 

35 000 (C XX 0.71 inter) 
40000 (C XX 0.71 hot) 

Porb 


8.3 d 


-'^orb 


34.9 Ro 



Notes. The different temperatures of the hotspot correspond to differ¬ 
ent simulations labelled (cold, inter, hot), ‘xx’ stands for ‘Id’ or ‘hd’ 
and designates different dust-to-gas ratios (see text). and Aorb are 
the orbital period and separation, respectively. Refer to Table 2 for fur¬ 
ther explanations. The hotspot radius is calculated using the black body 
relation (Sect. 2.4.3) and the estimate given by Eq. 3 is shown for com¬ 
parison in parenthesis. The density profile of model C is used for all 
three models (see Sect. 2.1). 


We can estimate the size of the hotspot (rhs) as being iden¬ 
tical to the size of the impacting stream (rstream), which can in 
turn be identified with the size of the stream at the Lagrangian 
point Xi. According to Pringle & Wade (1985), it can be ap¬ 
proximated by 


Lhs Atream 


1 Cs 

2 Gorb 


(3) 


where Cs is the sound speed of the material leaving the donor star 
at Xi and Qoib is the orbital angular velocity. 

The effective luminosity felt by the material escaping the sys¬ 
tem depends on its location along the spiral. In the cone with 
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opening 0flow above the hotspot (Fig. 1), the escaping material 
is irradiated by both the gainer and the hotspot^. But outside 
this cone, the hotspot luminosity does not contribute any more. 
We account for this effect by using an effective hotspot luminos¬ 
ity, corresponding to the actual value reduced by a geometrical 
factor y. This geometrical factor is obtained (Appendix B) by as¬ 
suming that the hotspot luminosity is diluted over a strip-shaped 
surface encircling the star (blue band in Fig. B.l). The effective 
flux from that surface reads 

Feff = Fg H- yFhs, (4) 

where Fg = L,/??!* and Fhs = Lhs/,^h.s are the surface fluxes of 
the gainer star and the hotspot, respectively (see Appendix B for 
details on the calculation of the factor y). With the quantities 
given in Table 1, we find y ~ 0.05 for model C. 

The use of a time-independent configuration, with an effec¬ 
tive flux that constantly irradiates all the surrounding material, 
is justified by the fact that the cooling time-scale (as computed 
by Cloudy; see Fig. 5) becomes longer than the orbital period at 
distances r « 3.8 x 10'^ cm (0.25 au) for model C. Due to the 
high outflow velocity (up to 900 km s"'; Fig. 3), the time needed 
for the flow to reach that distance is short (a; 13‘?93 = 1.67Porb 
at the end of the non-conservative phase). Therefore, the out¬ 
flow properties (like temperature, emission spectrum, ...) will 
not change too much over one orbital period, and several phases 
of irradiation (as the material passes in front of the hotspot) can 
be averaged into a constant irradiation. 

The spectrum of the gainer star is obtained by interpolating 
in the ‘TLUSTY’ grid of Lanz & Hubeny (2007) to match the 
effective temperature given in the first column of Table 1 . The 
donor-star spectrum is taken as a black-body spectrum (see up¬ 
per panel of Fig. 7). For all the components constituting the total 
incident spectrum (donor, gainer, hotspot), we use the metallic- 
ities listed in Table D.l. The donor has by definition the same 
abundances as the material ejected and, for the selected models, 
the gainer star is entirely covered by the donor’s material. 

2.4.2. Distance 

To relate absolute fluxes (from the models) with predicted ap¬ 
parent magnitudes for e.g., the infrared excesses, the distance to 
the system needs to be fixed. We choose a typical distance from 
the distribution of measured parallaxes for Algols. In Fig. 4, 
we present the distance distribution for 152 objects [join entries 
in the Budding et al. (2004) and Flipparcos (ESA 1997) cata¬ 
logues], which correspond to about 10% of all Algols in the 
Budding et al. (2004) catalogue. To evaluate which among these 
systems might be detectable by IRAS, Fig. 4 also presents the 
25 pm flux of the dusty simulation C_ld_0.71_inter, integrated 
over the whole nebula, as a function of distance (dashed green 
line). Only the infrared sources with a flux larger than the IRAS 
5cr detection limit at 25 pm (blue horizontal dotted line) can 
be detected, so that a system with the same IR emission as our 
model C_ld_0.71_inter must be within 400 pc of the Sun for its 
IR excess to be detected, and this includes most of the Algols 
with a known parallax. Based on this result, a typical distance 
of 300 pc is adopted to convert model fluxes into apparent fluxes 
and magnitudes, as will be listed in Sects. 3 and 4.2.2. To give 
an idea, (3 Lyrae is located at a distance of 294 pc while W Ser is 

^ We assume that the donor star does not contribute to the irradiation 
of the outflowing material, because its luminosity is much lower than 
that of the gainer (see Table 1). This, however, may not remain true in 
Algols with a fainter gainer of spectral type A. 
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Fig. 4. Distance distribution for observed Algols. The sample of Al¬ 
gols is from Budding et al. (2004). Parallaxes are from the Hipparcos 
catalogue (ESA 1997). The dashed green curve represents the flux at 
25 pm (right-hand scale) as a function of distance for the dusty simula¬ 
tion C_ld_0.71_inter. The flux is integrated over the whole nebula. The 
dotted blue line is the IRAS 5cr limit at 25 pm. The WISE 5cr limit for 
band 1V4 is 2600 pJy and is always below the dashed green curve. 


estimated to lie at 1428 pc (Hipparcos catalogue, van Leeuwen 
2007). Of course, the apparent fluxes of the model SEDs can 
be easily rescaled for comparison with a system at a different 
distance. 

2.4.3. The temperature of the hotspot 

To assess the dependence of our results on the SED of the 
hotspot, we use three different hotspot temperatures Ths in model 
C with /c = 0.71 (see Tables 1 and 2). 

In the first model (C_xx_0.71_cold), we assume that the 
hotspot emits with a B-star spectrum at a temperature of 
30000 K^. In the second model (C_ _xx_0.71_inter), we as¬ 
sume a temperature of 35 000 K and use an O-star spectrum 
(also interpolated from the TLUSTY grid). In the third model 
(C_xx_0.71_hot), we use an O-star spectrum with a temperature 
of 40000 K. These three temperatures are based on observed 
hotspot values in Algols with similar properties (stellar masses 
and orbital periods; van Rensbergen et al. 201 1), and the corre¬ 
sponding results are presented in Sect. 3.3. These three models 
allow to clearly disentangle the effect of the hotspot properties 
on the emergent SED. 

In order to satisfy the Stefan-Boltzmann law Lhs = 

(cf. Eq. B.2), the radius of the hotspot is treated as a free pa¬ 
rameter (Lhs and Ths are fixed), and compared to the estimate 
given by Eq. 3 (see Table 1). Although Eq. 3 can only give 
a rough approximation of the radius of the hotspot, we see a 
good match between this later value and the one imposed by the 
Stefan-Boltzmann law. 


2.5. Setup of Cloudy and Skirt 

We use version 13.02 of Cloudy, last described by Eerland et al. 
(2013), to calculate the ionisation of the gas close to the gainer 
star and the emission from the dust grains in the outflow at 

^ In our models, the hotspot is a region of the star’s atmosphere with 
a higher temperature and we therefore use a synthetic star spectrum 
instead of a black-body spectrum. 
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larger distances. Cloudy calculates the thermal and ionisation 
balance with a simultaneous treatment of radiative transfer. For 
the latter, an escape probability formalism is used. In our sim¬ 
ulations, the photon energy scale is discretised with a resolu¬ 
tion of d/A/i = 2 000 down to d = 0.152 nm'^. Cloudy cal¬ 
culates the strengths of ^ 10^ spectral lines. We found that 
the number of energy levels that are used for the Fe ii ion can 
have a strong impact on the resulting temperature profile and 
thus on the transition radius r^, and also on the emitted spec¬ 
trum. Therefore, we used the detailed Fe ii model that was devel¬ 
oped by Verner et al. (1999) and implemented in Cloudy. To de¬ 
fine the spatial step size. Cloudy uses adaptive logic. However, 
for the models C_df_0.7 l_hot, C_ld_0.7 l_hot, C_df_l .00_inter, 
and C_ld_1.00_inter, we had to define a lower limit for the step 
size, of the order of 10^ cm, to reach convergence. 

The ejection of matter triggered by the hotspot is assumed to 
occur preferentially in the orbital plane and is assumed to be seen 
edge-on. In our spherical ID simulations with Cloudy, we take 
that into account by adjusting the covering factor (equivalent to 
the factor /c introduced in Sect. 2.1, Eq. 1). In Appendix A we 
estimate the covering factor (or, equivalently, the opening angle) 
based on the assumption of isotropic radiation from the hotspot. 
To study the effect of the choice of the opening angle on the 
results, we performed simulations with different covering factors 
(the three digits in the label of each simulation: ‘x.xx’). 

To assess the systematic uncertainties due to the ID geom¬ 
etry and the simplified approach used by Cloudy for the radia¬ 
tive transfer, we use the Skirt code, which simulates the general 
3D continuum transfer for dust, including absorption, multiple 
anisotropic scattering, and thermal emission (Baes et al. 2003, 
2011; Camps & Baes 2015). The 3D dust-distribution is repre¬ 
sented by a disc that is centred on the gainer star and truncated at 
its inner edge at the dust condensation radius. The thickness of 
the disc increases linearly with distance as defined by the open¬ 
ing angle 0flow (see Appendix A), and the density inside the disc 
is assumed to be a function of distance only (see Fig. 3). The spa¬ 
tial grid consists of 200 concentric shells. The density is set to 
zero in the directions where the polar angle is outside the range 
n/2 - 0flow/2 < 0 < njl + 6flow/2, and the extent of the grid cells 
is adapted accordingly. The positions r, of the grid cells along 
the radial coordinate follow a power-law distribution of the form 

L/+1 = Ttr -L j _ ~ 

where N > i > 0 and ro - 0. ra is the dust formation radius 
(212 AU for the model C_ld_0.71_inter; Sect. 3.1), N = 200 
the number of grid points, and Tmax is chosen in such a way that 
the density drops below the ISM density at the outer edge of 
the grid. R = 800 is approximately the ratio of the sizes of the 
outermost and innermost grid cells, (r^r - rAr-i)/(r 2 - rj). The 
functional form of Eq. 5 is motivated by the steeply decreasing 
density profile. The ratio of the widths of two adjacent grid cells 
is held constant. This implies that the width of the grid cells in¬ 
creases from the centre towards the outer edge. Therefore, the 
high-density regions close to the centre are covered by a larger 
number of cells than the outer low-density regions. This allows 
us to save computing time. The first term on the right-hand side 
of Eq. 5 is introduced because there is a dust-free region around 
the centre, which does not have to be resolved in the Skirt simu¬ 
lations. The composition and size distribution of the dust are the 


/t/A/l = 333 between Cloudy’s high-energy limit at 12.4 fm and 
0.152 nm. 



Fig. 5. Temperature profile (solid red line) and cooling time-scale (dot¬ 
ted blue line; see text) of the material surrounding the star as a function 
of distance to the star (model C_df_0.71_inter). The dot-dashed cyan 
line represents the dust condensation temperature (1400 K). 

same as in the Cloudy models. In the radiative-transfer Monte- 
Carlo method used in Skirt, the radiation field is represented 
by discrete photon packages. We used 10® photon packages per 
wavelength in our simulations. Our wavelength grid consists of 
1 000 wavelength points with logarithmic spacing between 1 nm 
and 500 pm, which allows convergence within a reasonable com¬ 
putation time. 

The different parameters used in the simulations performed 
with Cloudy and Skirt are summarised in Table 2. 

3. Results 

3.1. Temperature and emissivity profiles 

The first quantity of interest is the outflow-temperature profile, 
because dust can only form if it falls below the condensation 
temperature. Eigure 5 shows the temperature profile for the 
model C_df_0.71_inter (see Table 2). The temperature drops 
below the dust condensation temperature at r^ = 3.2 x 10*® cm 
(ai 212 AU), which is small compared to the extent of the out¬ 
flow Tcut ~ 6.1 X 10*^ cm (see Sect. 2.1). Of course, this value 
depends on the covering factor and incident spectrum. 

The temperature profile shown in Eig. 5 is a decreasing func¬ 
tion of distance in the inner regions, with temperatures drop¬ 
ping from 17 000 K to » 180 K in the first » 1 x 10*^ cm 
(»;6 700 AU) from the star’s centre. Then, due to the steeply 
decreasing density, the temperature increases again up to about 
300 K. 

Two other interesting quantities to follow are the profiles 
of the volume emissivity at some specific wavelengths (espe¬ 
cially for the investigation of dust), like the central wavelengths 
for some widely-used IR-photometry systems (2MASS, WISE; 
Cutri et al. 2003, 2012), and the ionisation fraction profiles of 
some ions that are observed in the emission-line spectra of many 
Algol-related systems (Ee ii. Si iv, etc...). The ionisation fraction 
for an ion X z is defined as luxz = nxz/nx,total, where ux.totai is the 
total number density of element X. These profiles tell us where 
most of the flux comes from. 

The top panel of Eig. 6 shows the emissivity profiles of the 
2MASS Ks, WISE W2 and W4 bands for the dusty simulation 
with fc = 0.71 and O = 2 x 10 ® (model C_ld_0.71_inter). Most 
of the flux comes from the inner part of the system, but the in- 
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Table 2. Parameters of the simulations performed with the Cloudy and Skirt codes. 


code 

Cloudy 

Skirt 

fc or ^flow 

0.60/0.71/1.00 

o 

o 

density 

plateau 

normal 

normal 

chemistry 

gas 

gas -H dust 

gas 

gas + dust 

dust 

O 

— 

2 X 10-5 

— 

2 X 10-5 

4 X 10-5 

2 X 10-5 

4 X 10-5 

label 

C_df_p_0.71 


A_df_0.71 

B df 0.71 

C df 0.71 cold 
C_df_x.xx_inter 
C_df_0.71_hot 

A_ld_0.71 

B Id 0.71 

C Id 0.71 cold 
C_ld_x.xx_inter 
C_ld_0.71_hot 

C hd 0.71 

C SK Id 

C SK hd 


Notes, ‘plateau’ (models labeled ‘p’) indicates that the simulation has a constant density in the inner part (see text, Sect. 2.2). O is the dust-to-gas 
ratio, /c the covering factor (used for the Cloudy simulations) and 0flow the opening angle (used in the Skirt simulations to define a 3D dust 
distribution). In all models, the hotspot temperature is set to Th. = 35 000 K (i. e. ‘inter’), except for the models with /(, = 0.71, for which we 
applied three different hotspot temperatures (‘cold’: Ths = 30000 K; ‘inter’: Th., = 35 000 K; ‘hot’: Ths = 40000 K) with different spectra (see 
Sect. 2.4.3). ‘A’, ‘B’ and ‘C’ refer to the evolutionary phase (see Fig. 2 and Table 1), ‘x.xx’ is a three-digit number for the covering factor, ‘df’ 
stands for dust-free, ‘Id’ for low dust (i. e., O = 2 x 10“^), ‘hd’ for high dust (i. e., cD = 4 x 10“'^). 


frared excess in these specific bands comes from the outer dusty 
regions at r ^ 2 x 10*^ cm. We indeed clearly see that W4 
dominates W2 and 2MASS in this outer region, which is not 
the case at smaller radii. In the outer region, the emissivity is 
low (of the order of 10“^^ ergcm^^ s“') but the total emission, 
which is the integral over the volume, remains large. The in¬ 
frared excess in the W4 band is therefore visible as a ‘bump’ in 
the infrared part of the SED from ~10 pm towards longer wave¬ 
lengths (Fig. 9). Due to the rapidly decreasing density profile, 
only the dust the closest to the star contributes to the infrared ex¬ 
cess, especially at the lowest wavelengths. For this simulation, 
as well as for the dust-free one, the inner part is totally ionised 
(Hit, Feiii and higher). Although Feii is mainly present in the 
outer part of the system, the formation of Fe ii lines may par¬ 
tially happen close to the star and present some line variability 
there (see Sect. 3.3 and bottom panel of Fig. 6). In contrast, the 
contribution from Fe ii in the outer parts comes from a smooth, 
roughly time-steady density distribution and is probably not vari¬ 
able. The resulting observed line will be the sum of the two con¬ 
tributions (variable inner part, steady outer part). 

3.2. Synthetic SEDs 

In this section, we inspect the output spectra in the different 
regimes, from ultraviolet to infrared, and the impact of the dif¬ 
ferent parameters. 

3.2.1. Optical-UV 

The top panel of Fig. 7 shows the emergent spectral energy 
distributions (SEDs) for the dust-free models C_df_0.60_inter, 


C_df_0.7l_inter, C_df_0.7l_cold and C_df_l .00_inter (the cor¬ 
responding SEDs for models that contain dust are nearly identi¬ 
cal in this spectral region and are therefore not shown). They all 
show pronounced Balmer and Paschen recombination continua. 
The Balmer jump in emission leads to observable diagnostics, as 
further discussed in Sect. 4.3. 

In the SEDs of the models C_df_0.71, we notice a slight ab¬ 
sorption in the UV up to a wavelength of ~ 0.2 pm (also seen 
in Fig. 8 for the same model), which is mainly attributable to 
photoionisation close to the star, where the density is high. This 
has two consequences; i) this absorption will quickly disappear 
when the non-conservative phase will stop, because the outflow¬ 
ing material will rapidly dissipate (see Sect. 2.1); ii) the UV ab¬ 
sorption should show a phase dependence owing to the spiral 
structure of the outflow (this prediction is compared to observa¬ 
tions of W Ser in Appendix C). 

The bottom panel of Fig. 7 presents the SED for the differ¬ 
ent models (A, B, C) with a covering factor /c = 0.71 (cor¬ 
responding to fiflow = 90°) and Ths = 35 000 K. Model C ex¬ 
hibits the strongest Balmer and Paschen recombination continua 
and model A the most prominent UV absorption. As discussed 
in Sect. 4.3, The presence of the recombination continuum and 
emission lines, and the absence of stellar absorption lines are not 
typical of Algols, but are a common feature of symbiotic systems 
or planetary nebulae (Kenyon 1986). 

3.2.2. Far UV-continuum orbital modulation 

The comparison of the results from the model C_df_0.71_inter 
and the corresponding ‘plateau’ model C_df_p_0.71 (the lat¬ 
ter having a lower density in the inner region around the sys¬ 
tem) allows us to evaluate the impact of the azimuthal depen- 
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C_dlu0.71_inter 
C_dlU1.00_inter 
C_dfi0.71_cold 
C df iO.60 inter 



0.6 0.7 0.8 0.9 1.0 


Fig. 7. Top: Synthetic SEDs of the dust-free Cloudy models in the UV-optical regime. The input spectrum is the sum of the spectra of the 
gainer star, donor star, and hot spot with a temperature of 35 000 K (the latter diluted according to the prescriptions of Appendix B). The fluxes are 
computed for a system distance of 300 pc. Bottom: Same as top panel but for the different models A, B and C (at different evolutionary times). 
The input spectra are in black from bottom to top (A, B, C). 


dence of the outflow density, as the line of sight probes different 
parts of the outflow spiral structure. As discussed in Sect. 2.2, 
the ’plateau’ model mimicks situations corresponding to orbital 
phases where the line of sight does not enter the cone of opening 
0flow (Fig- 1), thus when the density close to the star is low. Fig¬ 
ure 8, which presents the SED in the UV regime for these two 
models, reveals that the ’plateau’ model has no far-UV absorp¬ 
tion, whereas the non-plateau models have. Thus, one expects 
an orbital modulation of the far-UV continuum as a consequence 
of the azimuthal distribution of the outflowing material close to 
the star. A confrontation with far-UV spectro-photometric ob¬ 
servations could thus possibly provide insights into the spatial 
structure of the outflow. 


3.2.3. Emission lines from weakly ionised species 

Figure 7 reveals the presence of many Fe n emission lines for the 
different models considered, which are best seen in the range 0.2 
- 0.3 pm but are also present in the optical. Figure 6 reveals that 
these lines mainly form in the outer part of the outflow, outside 
the Roche lobe. Detecting these lines may help to prove that 
some material is indeed escaping from the binary system. The 


absence of such lines in spectra of genuine Algols may thus be 
considered as an argument against the presence of outflowing 
matter in these systems. 

3.2.4. Infrared 

Figure 9 shows the characteristic silicate feature that is produced 
by the dust grains in the outer part of the outflow in our models. 
All models produce infrared excesses above the 5cr IRAS de¬ 
tection limit at 12 and 25 pm, but only the simulation with the 
largest dust-to-gas ratio has a 60 pm flux detectable by IRAS (for 
systems located as usual at a distance of 300 pc). In the ‘plateau’ 
simulation (C_ld_p_0.71, not shown in the figure), there is no 
dust emission, because the temperature falls below the conden¬ 
sation temperature at a location where the density is too low to 
produce any detectable infrared excesses. 

3.2.5. Synthetic photometry 

To ease the comparison of our models with observations. Ta¬ 
ble 3 presents, for all our models, the magnitudes in the Strom- 
gren uvby system and the corresponding colour indices, using the 
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Table 5. Magnitudes in 2MASS and WISE bands assuming the models are located 300 pc away from the sun. 


Model 


J 

H 


W1 

W2 

W3 

W4 


FwaIF-wi 

Fwi / Fwi 

A 

df 

0.71 


5.24 

5.46 

5.23 

5.31 

5.23 

5.01 

4.70 

2.39 

0.05 

0.14 

A_ 

Jd_ 

_0.71 


5.23 

5.46 

5.23 

5.31 

5.23 

4.19 

3.22 

2.41 

0.19 

0.29 

B_ 

_df 

_0.71 


5.52 

5.72 

5.45 

5.46 

5.27 

4.95 

4.58 

2.27 

0.06 

0.16 

B 

Id 

0.71 


5.52 

5.72 

5.44 

5.46 

5.28 

4.24 

3.35 

2.26 

0.19 

0.31 

C 

_df 

T3 

1 

o 

5.67 

5.60 

5.51 

5.43 

5.35 

4.93 

4.58 

2.09 

0.06 

0.16 

C 

_df 

1.00 

jnter 

5.44 

5.68 

5.59 

5.77 

5.71 

5.46 

5.08 

2.80 

0.05 

0.14 

C 

Id 

1.00 

inter 

5.44 

5.68 

5.56 

5.77 

5.70 

4.61 

3.41 

2.80 

0.24 

0.30 

C 

_df 

0.71 

jnter 

5.39 

5.63 

5.44 

5.51 

5.41 

5.11 

4.71 

2.54 

0.06 

0.15 

C 

_df 

0.71 

_cold 

5.45 

5.67 

5.42 

5.45 

5.31 

4.99 

4.64 

2.37 

0.06 

0.16 

C 

Id 

0.71 

inter 

5.39 

5.63 

5.44 

5.51 

5.40 

4.49 

3.27 

2.54 

0.21 

0.26 

C 

Jd_ 

0.71 

_cold 

5.44 

5.67 

5.43 

5.46 

5.32 

4.21 

3.20 

2.39 

0.22 

0.33 

C 

M 

0.71 

_inter 

5.39 

5.64 

5.44 

5.49 

5.17 

-0.21 

-1.99 

2.54 

26.43 

19.43 

C 

Id 

0.71 

hot 

5.37 

5.63 

5.45 

5.53 

5.41 

4.35 

3.17 

2.59 

0.24 

0.30 

C 

_df 

0.71 

_hot 

5.37 

5.63 

5.45 

5.53 

5.42 

5.14 

4.74 

2.60 

0.06 

0.15 

C 

Jd_ 

0.60 

jnter 

5.38 

5.59 

5.31 

5.34 

5.21 

4.20 

3.12 

2.27 

0.21 

0.29 

C 

_df 

0.60 

jnter 

5.38 

5.59 

5.31 

5.34 

5.23 

4.93 

4.60 

2.28 

0.05 

0.15 

Stellar SED (C) 

5.76 

5.77 

5.76 

5.76 

5.75 

5.73 

5.74 

2.41 

0.04 

0.11 


Notes. FjIFk^, Fy/^IFyfi and are the flux ratios. 


Table 3. Magnitudes in Stromgren (uvby) bands and corresponding 
colour indices. 


Model 


u 

V 

b 

y 

b-y 

m\ 

Ci 

A_ 

_df 

_0.71 


5.08 

5.31 

5.22 

5.20 

0.02 

0.07 

-0.32 

A_ 

Jd_ 

_0.71 


5.09 

5.31 

5.22 

5.20 

0.02 

0.07 

-0.32 

B 

df 

0.71 


5.33 

5.64 

5.54 

5.49 

0.05 

0.05 

-0.41 

B 

Jd_ 

0.71 


5.33 

5.64 

5.54 

5.49 

0.05 

0.05 

-0.41 

C 

df 

1 

o 

5.75 

5.68 

5.67 

5.71 

-0.03 

0.04 

0.07 

C 

_df 

1.00 

jnter 

5.21 

5.48 

5.42 

5.34 

0.08 

-0.02 

-0.33 

C 

Jd_ 

1.00 

jnter 

5.21 

5.48 

5.42 

5.34 

0.08 

-0.02 

-0.33 

C 

df 

0.71 

inter 

5.08 

5.37 

5.27 

5.15 

0.11 

-0.02 

-0.39 

C 

_df 

0.71 

_cold 

5.17 

5.48 

5.39 

5.35 

0.04 

0.05 

-0.40 

C 

Jd_ 

0.71 

jnter 

5.08 

5.37 

5.27 

5.15 

0.12 

-0.01 

-0.39 

C 

Id 

0.71 

_cold 

5.18 

5.48 

5.39 

5.35 

0.04 

0.05 

-0.39 

C 

Jid 

0.71 

Jnter 

5.08 

5.37 

5.27 

5.16 

0.12 

-0.01 

-0.39 

C 

Jd_ 

0.71 

Jot 

5.07 

5.34 

5.25 

5.14 

0.12 

-0.03 

-0.36 

C 

df 

0.71 

Jot 

5.07 

5.34 

5.25 

5.14 

0.12 

-0.03 

-0.36 

C 

Jd_ 

0.60 

jnter 

5.15 

5.42 

5.33 

5.33 

-0.01 

0.10 

-0.35 

C 

df 

0.60 

jnter 

5.15 

5.42 

5.33 

5.33 

-0.01 

0.10 

-0.35 

Stellar SED (C) 

5.38 

5.28 

5.28 

5.33 

-0.05 

0.05 

0.10 


Notes, nil = (v- b) - (b - y); Ci = (u- v) — (v- b) 


passbands and zero points of Bessell (2011) and assuming that 
the models are located 300 pc away from the sun. As expected, 
the systems with a strong Balmer continuum emission appear 
with a negative value of ci - {u-v)-{v- b). Figure 10 presents 
our models in the {m\ Q ; ci,o] diagram compared with data from 
the Long-Term Photometry of Variables project (Manfroid et al. 
1991). Only comparison stars (labelled A, B, or C in their cata¬ 
logue) have been retained to define the main sequence. The ci,o 
and OTi 0 indices are not affected by interstellar reddening and are 


Table 4. Magnitudes in SDSS bands and corresponding colour indices. 


Model 


u' 

9’ 

r' 

i' 

z' 

u' - g' 

g'-r' 

A df 0.71 


4.83 

5.06 

5.23 

5.40 

5.86 

-0.22 

-0.17 

A Id 0.71 


4.84 

5.06 

5.23 

5.40 

5.86 

-0.22 

-0.17 

B df 0.71 


5.09 

5.36 

5.51 

5.68 

6.15 

-0.27 

-0.15 

B Id 0.71 


5.09 

5.36 

5.51 

5.68 

6.15 

-0.27 

-0.15 

C_df_p_0.71 

5.43 

5.51 

5.81 

6.04 

6.23 

-0.08 

-0.30 

C df 1.00 

jnter 

4.96 

5.22 

5.25 

5.44 

5.71 

-0.26 

-0.03 

C Id 1.00 

jnter 

4.96 

5.22 

5.25 

5.44 

5.71 

-0.26 

-0.03 

C df 0.71 

jnter 

4.84 

5.08 

5.24 

5.44 

5.89 

-0.25 

-0.16 

C df 0.71 

_cold 

4.93 

5.22 

5.39 

5.57 

6.02 

-0.28 

-0.18 

C Id 0.71 

jnter 

4.84 

5.08 

5.24 

5.44 

5.89 

-0.25 

-0.16 

C Id 0.71 

_cold 

4.94 

5.22 

5.39 

5.57 

6.02 

-0.28 

-0.17 

C hd 0.71 

inter 

4.84 

5.08 

5.24 

5.44 

5.89 

-0.25 

-0.16 

C Id 0.71 

Jot 

4.82 

5.05 

5.15 

5.42 

5.87 

-0.23 

-0.10 

C df 0.71 

Jot 

4.82 

5.05 

5.15 

5.42 

5.87 

-0.23 

-0.10 

C Id 0.60 

jnter 

4.90 

5.17 

5.36 

5.55 

5.99 

-0.27 

-0.19 

C df 0.60 

jnter 

4.90 

5.17 

5.36 

5.55 

6.00 

-0.27 

-0.19 

Stellar SED (C) 

5.30 

5.33 

5.68 

5.97 

6.21 

-0.03 

-0.35 


defined by (Stromgren 1966): 

ci.o =ci - 0.2{b - j/) = [(« b)] - 0.2(b - y) (6) 

mi,o =mi + O.mb -y)^ [(n - b) - (b - y)] + O.mb - y). (7) 

While our input spectrum falls within the main-sequence pop¬ 
ulation (Stromgren 1966), our models considerably deviate to¬ 
wards lower ci_o- Deviations in oti q are also observed for some 
models. In addition to the Stromgren photometry. Table 4 shows 
the SDSS magnitudes and two color indices, which we calcu¬ 
lated using the passbands provided by Moro & Munari (2000) 
and the zero points from Fukugita et al. (1996). We also present, 
in Table 5, the 2MASS and WISE band magnitudes that will be 
used in Sect. 4.2.2. WISE magnitudes are computed according 
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Fig. 6. Top panel: Emissivity profiles (solid lines) and relative ion¬ 
isation (dashed lines) for the simulation C_ld_0.71_inter (blue solid 
line: 2MASS K^; green and red solid: WISE IV2 and W4, respectively; 
green dashed: Hit fraction; blue dashed: Feii fraction). The vertical 
cyan line represents the Roche lobe radius for the gainer star. Bottom 
panel: Si iv (solid lines) and Fe ii fraction (dashed lines) for the high- 
temperature hotspot simulations (blue: model C_df_0.71_inter; red: 
C_df_0.71_hot). 
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Fig. 9. Same as Fig. 7 but for the IR region, including the Skirt 
model C_SK_ld. For clarity, spectral lines have been removed from 
the Cloudy SEDs. The red arrows are 5 <t IRAS upper limits (in the 
faint-source catalogue), namely 0.2 Jy for the 12 and 25 |.tm bands and 
0.28 Jy for the 60 tim band. 
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Fig. 10. Stromgren photometry [mlo ; cIq] diagram for our mod¬ 
els (crosses) and for comparison normal main-sequence stars (red dots) 
from the Long-Term Photometry of Variables project (Manfroid et al. 
1991, to be discussed in Sect. 4.3). From the latter, only data obtained 
with the Danish telescope (system 7) were retained. For clarity, the 
models are not labelled and only their general properties are indicated 
in blue. We refer to Table 3 for the values. We added observations of 
W Ser, B[e] and Be stars for comparison (green dots). The blue dots are 
the average values for the latter objects when several observations are 
available. 



according to ADPS® (See also the “Explanatory Supplement to 
the 2MASS All Sky Data Release and Extended Mission Prod¬ 
ucts”^). 


Fig. 8. Synthetic SEDs of dust-free models in the UV regime. The 
solid black line in the upper panel represents the input stellar spec¬ 
trum. The lower panel shows the difference between the output and 
input spectra for the two models identified by the line labels. 

to Tokunaga & Vacca (2005), based on transmission curves from 
the Explanatory Supplement to the WISE All-Sky Data Release 
Products (Cutri et al. 2012)^. 2MASS response curves are set up 

^ http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/index. html, 2014-02-27, Sect. IV.4.h, Figure 5a and Table 3 


3.2.6. Radio emission: H i-21 cm and rotational CO lines 

The 21 cm line of neutral hydrogen (Hi-21 cm) and molecular 
rotational lines such as COi_o d = 2.6806 mm and CO 2-1 d = 
1.3004 mm may help to map the cold (neutral) outer regions of 
the outflow. Although the 21 cm line may be completely blended 

^ http://ulisse.pd.astro.it/Astro/ADPS/Systems/ 
Sys_new_021/index_n021 .html, 2014-03-05 
’ http://WWW.ipac.caltech.edu/2mass/releases/ 
allsky/doc/explsup.html, 2014-03-05, Sect. VI.4.a 
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^C) 


Fig. 11. [Siiv] 1394 and 1403A doublet for three different configura¬ 
tions. 


with ISM lines, the high speed of the ejected material (with ter¬ 
minal velocities up to «i900 km s“^) may help to distinguish 
between the circumstellar and ISM contributions. In our mod¬ 
els, the CO abundance is very small and does not lead to any 
observable lines. The Hi-21 cm line is also very weak, possibly 
because of the low densities far from the star, in the neutral hy¬ 
drogen region. The Hi-21 cm line luminosity does not exceed 
a few 10“*^ erg s^' for both the dust-free and dusty simulations, 
corresponding to a few pjy at a resolving power of 2 000 and 
a distance of 300 pc. The outer part of the outflow is therefore 
probably not detectable in these radio bands. 

3.3. Models with strong ionisation 

In this section, we investigate the impact of the hotspot temper¬ 
ature and SED on the emergent SED as well as on the formation 
of dust and emission lines. 

A hotspot SED corresponding to a higher temperature im¬ 
plies a higher number of ionising photons and thus larger abun¬ 
dances of species such as Siiv, Nv or Civ. Eigure 11 shows 
the [Siiv] doublet at 139.4 nm and 140.3 nm. There is no dif¬ 
ference between the incident and emerging SEDs for the model 
with Ths = 30000 K (’C_df_0.71_low’). With a hotter hotspot 
and star, emission lines from strongly ionised elements can form 
(although still only barely noticeable for Ths = 35 000 K; model 
’C_df_0.71_inter’). 

The bottom panel of Eig. 6 shows the relative abun¬ 
dance of Siiv and Eeii for the models C_df_0.71_inter and 
C_df_0.71_hot. We find that Si iv is restricted to the Roche vol¬ 
ume around the gainer star. We therefore expect a strong line- 
profile variability due to the orbital motion, as well as eclipses 
by the donor and the gainer stars as observed in some W Ser sys¬ 
tems (see Appendix C). Unfortunately, we cannot simulate these 
effects in our models, which do not resolve the inner binary sys¬ 
tem. Due to the short cooling time-scale close to the star (Eig. 5), 
only the material above the hotspot (within the cone of opening 
angle 0flow in Eig. 1) may contain highly ionised species. The low 
hotspot-temperature models produce much less Si iv, but more 
Ee II than the high-temperature ones. 

The observed presence of [Siiv] emission in the Algol- 
related systems W Ser (Sect. 4.1 and Appendix C) does not nec¬ 
essarily imply systemic mass loss however, because it probes 
the inner part of the nebula, which may still be gravitation- 
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Fig. 12. Skirt intensity map at 12 pm (top), 25 pm (middle) and 60 pm 
(bottom) for O = 2 X 10“^ (model C_SK_ld). To ease the comparison, 
the same limits are used for the three bands for the ‘x’ and ‘y’ axes. The 
angular size and intensity are given for a system at 300 pc. The figure 
is zoomed on a box 12 000 x 8 516 AU^ wide. 


ally bound. However, if present, the variability of lines from 
such highly ionised species may help to understand the ejection 
mechanism (Guinan 1989). Indeed, since the ionised elements 
form close to the hot region (the ejection point in the hotspot 
paradigm), we should be able to constrain the position of the 
UV source and check whether it is compatible with the location 
of the accretion point. On top of that, radial-velocity measure¬ 
ments may constrain the global motion of this material (expelled 
from the accretion region or revolving around the accreting star, 
...). 

With increasing hotspot temperatures, only the outermost re¬ 
gions contain dust. Thus the following dichotomy is expected: 
[dust/weak ionisation] or [no dust/strongly ionised elements], 
which is supported by observations as discussed in Sect. 4.2. 

3.4. Parameter study: geometry - Skirt simulations - 
dust-to-gas ratio 

In this section, we compare ID results from Cloudy with re¬ 
sults from the 3D dust radiative-transfer code Skirt. When the 
codes are applied to the same environment (as described by the 
stellar parameters, the opening angle of the stream, ...), a very 
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good agreement between the ID and 3D results is found (see 
Fig. 9: compare the dashed green line and dotted long-dashed 
cyan line). This agreement validates the results obtained with 
Cloudy and described so far. 

As described in Sect. 2.5, we calculated models with differ¬ 
ent covering factors /c representing different outflow geometries. 
Figure 7 shows the impact of the covering factor for three con¬ 
figurations: the spherical emission case (/c = 1) and two cases of 
a disc-like emission (/c = 0.71 and /c = 0.6). For each covering 
factor, we recompute the density profile according to Eq. 1 so 
that the total mass in the outflow is conserved. The UV/optical 
absorption increases with decreasing covering factor, because 
the densities increase, but the effect is modest. With the higher 
covering factors, the IR flux is significantly reduced because the 
dust condensation radius being larger, less dust contributes to the 
infrared emission (Fig. 9). 

Figure 12 shows a 2D map of the edge-on system at three 
different wavelengths (12, 25, and 60 pm) produced with Skirt 
using the same physical properties as in model C_ld_0.71_inter 
(same opening angle and dust-to-gas ratio <t) = 2 x 10“^; 
model C_SK_ld). We restrain our investigation to regions 
where the flux is of the order of actual detection limits, set at 
0.1 mJy arcsec^^ («; 1/10 of current observational thresholds in 
the mid- to near-lR for the VLT/VISIR instrument; Lagage et al. 
2004). We do not expect a strong contribution from the cosmic 
infrared background (CIB) at the selected wavelengths. Schlegel 
et al. (1998) show that at 100 pm (their lowest wavelength), the 
CIB value is only =s3.3 mJy arcsec^^ in the galactic plane and 
much less closer to the galactic pole and at shorter wavelengths. 
The speckle effect (coarseness observed especially far from the 
centre of the box) seen in Fig. 12 is due to the smaller number 
of photon packages sent from the outer regions compared to the 
inner regions; it is thus numerical noise. The edges of the disc¬ 
like outflow are too cold to be seen at 12 pm and only the longer 
wavelengths can probe these outer regions. However, the flux 
decreases at longer wavelengths and drops below the adopted 
detection limit. The ‘X’ shape is due to the opening angle of 
the disc-like outflow, seen through the optically thin outer part. 
We remind that, in our models, there is no material outside the 
opening angle 6>flow, whereas in a real system we expect a much 
smoother transition, with the highest density on the binary plane 
and a lower density at higher latitudes. Therefore, the ‘X’ will 
be less marked. The inner hot region, where dust cannot form, 
has a diameter of a; 425 AU corresponding to » 1'.'4 (18 pixels) 
in Fig. 12. For comparison, the binary system has a separation 
of less than 1 AU and is thus far from being resolved. 

For the models (C_hd_0.71 and C_SK_hd) with a higher 
dust-to-gas ratio (O = 4 x 10“^, typical of the ISM) and a cover¬ 
ing factor fc = 0.71, the 12 pm peak rises up to 60 Jy in both 
the Cloudy and the Skirt simulations, as compared to 0.6 Jy for 
the model C_ld_0.71_inter displayed in Figs. 9 and 12. A ratio 
of 100 between the 12 pm fluxes of models differing by a fac¬ 
tor of 200 in their dust-to-gas ratios may seem surprising when 
the dust shell is optically thin, as is the case here. However, one 
should first subtract the stellar contribution from the observed 
12 pm flux of model C_ld_0.71_inter, and then take the ratio 
with the flux from model C_hd_0.71. The expected ratio of 200 
is found then. Clearly, the uncertainty on the dust-to-gas ratio in 
the Algol environment has a strong impact on the present predic¬ 
tions. 


4. Confrontation to observations 

4.1. Stellar samples 

In this section, we present the samples of objects investigated for 
the presence of systemic mass loss as predicted by the models. 
The observable diagnostics extracted from the models and dis¬ 
cussed in Sect. 3 are not typical of Algol systems, but share sim¬ 
ilarities with the properties of Algol-related systems like W Ser- 
pentis and p Lyrae stars, B[e] and Be stars. Symbiotic Algols and 
Double Periodic Variables. 

Algols. The definition of this class adopted by the Simbad 
database^ In the present study, we call “genuine Algols” systems 
that are detached or semi-detached, showing the Algol paradox 
in the absence of UV emission lines (that are used as a property 
of W Ser systems as discussed next). We created our own sam¬ 
ple of Algols based on the catalogues of Brancewicz & Dworak 
(1980) containing 701 Algols® and (Budding et al. 2004) con¬ 
taining 435 Algols^®. In these two catalogues, the classification 
‘Algol’ is based on the individual mass ratio. 

W Serpentls and /? Lyrs are the prototypes of two subsets 
of Class-Il'^ Algol systems. The accreting star in these objects 
may be obscured by a thick accretion disc (see the review of 
Guinan 1989). They correspond to rapid mass-transfer systems 
most likely during case B, or BB for p Lyr itself (second phase 
of mass transfer following the case-B mass-transfer phase). The 
main property of W Ser systems is the presence of strong ultravi¬ 
olet emission lines ([Civ], [Cii], [Hen], [N v], [Siiv]...; see Ap¬ 
pendix C) seen at every phase (Plavec 1980, 1982, 1988, 1989; 
Guinan 1989). Our sample of W Ser is built on the catalogue of 
31 peculiar emission-line Algols (PELAs; Gudel & Elias 1996) 
after removing from that list those systems with no evidence for 
the presence of these UV emission lines. The strong ultraviolet 
emission lines from highly ionised species observed in W Ser 
systems may be identified with those appearing in our model 
with the highest-temperature hotspot (C_df_0.7l_hot). However 
these spectral signatures cannot confirm the presence of material 
escaping the system since they form inside the Roche lobe of 
the star (see Sect. 3.3). Moreover, in W Ser systems, there are 
two contributions to the strongly ionised emission lines (Wei- 
land et al. 1995, see Appendix C): the boundary layer between 
the disc and the star, responsible for the broad emission seen at 
every phase and the hotspot, responsible for the highly phase- 
dependent narrow emission. There have also been reports of 
infrared emission that we will further investigate. Taranova & 
Shenavrin (1997) report infrared excesses at 3.5 and 5 pm possi¬ 
bly attributable to dusty material. Similarly, RY Set possesses an 
extended dusty nebula (2 000 AU wide; Grundstrom et al. 2007). 

Be stars are rapidly rotating (v>300 km s '), non-supergiant 
B stars showing strong Balmer emission lines or evidence of an 
ionised shell (see Rivinius et al. 2013 for a recent review). These 


* http: //simbad. u- strasbg. fr/simbad/ is “detached eclipsing 
binaries showing the Algol paradox (the more massive component of 
the system is the less evolved)”. However, many among the systems 
considered as Algols in the literature are actually semi-detached sys¬ 
tems. 

® http://edsare.u-strasbg.fr/viz-bin/Cat?II/150A 
http://edsare.u-strasbg.fr/viz-bin/Cat?J/A%2bA/417/263 
rapid mass transfer systems with mass ratio reversed (Deschamps 
et al. 2013) 
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stars form a dust-free, outwardly diffusing, disc. In single Be 
stars, the formation mechanism of these decretion discs is related 
to the rapid rotation of the star, but some Be stars are interacting 
binaries, likely Algols. In this case, the formation of the disc 
and the star’s rapid rotation could result from the RLOF accre¬ 
tion stream (Paczynski 1991; Colpi et al. 1991; Deschamps et al. 
2013). Although there is no catalogue of truly “Algol-like” Be 
stars yet available, Harmanec (2001) has compiled an extensive 
list of Be stars in binary systems that may serve as an interesting 
starting point. In that list, we identified five Algol-like systems 
of particular interest as it will be apparent from our subsequent 
analysis: <p Per (Gies et al. 1998), SX Aur (Linnell et al. 1988), 
HR 2142 (Peters 1983; Waters et al. 1991), n Aqr (Bjorkman 
et al. 2002) and CX Dra (Berdyugin & Piirola 2002). 


B[e] stars are B stars showing strong forbidden emission lines, 
Balmer emission lines and infrared excesses due to hot cir- 
cumstellar dust (Lamers et al. 1998). However, this class con¬ 
tains several poorly understood systems and includes young and 
evolved, single and binary stars. Therefore, different mecha¬ 
nisms at the origin of the gas surrounding B[e] stars may be 
involved, among which some may be related to mass transfer 
across the B[e] binary systems (see, e. g., Dunstall et al. 2012). 
There is no specific catalogue of binary B[e] stars. 


Symbiotic Algols. Although most symbiotic stars are assumed 
to host an accreting white dwarf (e. g., Mikolajewska 2007), it is 
not necessarily the case for all of them. For example, SS Lep 
(Blind et al. 2011) is a non-eclipsing symbiotic, showing the 
Algol paradox, composed of a red-giant donor and a main- 
sequence accreting star. This system shows the presence of a 
dusty circumbinary disc or envelope (Jura et al. 2001). The main 
difference with standard Algols lies in the mass transfer scheme, 
which is believed to be Wind Roche-lobe overflow (WRLOF; 
Mohamed & Podsiadlowski 2007) instead of standard RLOF. 
Due to the long period of the system (260^3), mass accretion 
onto the gainer is supposed to occur via an accretion disc. Un¬ 
fortunately, the properties of the hotspot on the edge of an ac¬ 
cretion disc are badly constrained, which makes it difficult to 
predict systemic mass-loss rates and to compare with observa¬ 
tional data. Moreover, it is not clear whether the hotspot mecha¬ 
nism can work with WRLOF, because of the lower mass-transfer 
rates. 


Double periodic variables (DPVs) are interacting binaries ex¬ 
hibiting two closely linked periodic variations, the shortest one 
being the orbital period (Mennickent et al. 2003, 2005, 2008), 
and the second one of unknown nature but probably linked to 
the variations in the strength of the disc wind (Mennickent et al. 
2012). They are considered as one specific evolutionary step for 
the more massive Algols. An interesting property of these ob¬ 
jects is the surprising constancy of their orbital period, which is 
not expected in Algols undergoing RLOF mass transfer with a 
mass ratio different from 1 (Garrido et al. 2013). An explana¬ 
tion could be that systemic mass loss compensates for the pe¬ 
riod increase caused by mass transfer from a less massive star 
to its more massive companion. Mennickent et al. (2008) al¬ 
ready emphasised the need for further investigation of these sys¬ 
tems. Some DPVs (the Galactic ones) already appear in our Al¬ 
gol or W Ser catalogues but most of these objects have been 
found thanks to surveys in the Small Magellanic Cloud and are 


therefore too far away to allow any clear detection of an infrared 
excess. We therefore do not include them in this study. 

4.2. Constraint on the dust formation: IR diagnostics 

In this section, we search for systems that exhibit infrared 
colours typical of dust emission. We select Class II and III ob¬ 
jects (flaring and quiescent Algols, respectively, with a mass ra¬ 
tio already reversed, according to the nomenclature of Paper I), 
for which we expect most of the mass to have been ejected, and 
dust to have formed. Class III objects should no longer show any 
signature of systemic mass loss, though, because of the short¬ 
living outflow, but they were nevertheless considered for com¬ 
pleteness. We extracted the IR data from IRAS, 2MASS and 
WISE catalogues. The use of the IRAS catalogues is motivated 
by their completeness in the 12 to 60 pm range. On the other 
hand, the choice of the WISE and 2MASS databases has been 
motivated by their better detection limits as compared to IRAS 
(IRAS sensitivities: 0.5 Jy at 12, 25 and 60 pm; WISE 5cr sen¬ 
sitivities: 120, 160, 650 and 2600 pJy at 3.3, 4.7, 12 and 23 pm, 
respectively; 2MASS lOcr at 7 < 15.0 mag, H < 14.2 mag, 
Xs < 13.5 mag). 

4.2.1. IRAS data 

We selected the data from the IRAS (InfraRed All-sky Survey) 
point'^ (Helou & Walker 1988) and faint*^ (Moshir & et al. 
1990) source catalogues. We retrieve 241 Algol entries with 
a matching IRAS source within a 2' positional tolerance. We 
search for IR excesses using the colour indices 

[/]-[;] = -2.5 log (8) 

where i and j correspond to the IRAS bands at 12, 25, or 60 pm, 
and Fjfl and E; o are the zero-magnitude fluxes. Rejecting spu¬ 
rious detections (contaminated by a nearby infrared source), 
we did not And any systems presenting a deviation larger than 
0.5 mag above the black-body line in the ([12]-[25],[12]-[60]) 
plane. 

4.2.2. WISE data 

2MASS'"^ and WISE'^ data are provided by Cutri et al. 
(2003) and Cutri et al. (2012), respectively. Eigure 13 
shows the flux ratios [WISE Fw 4 (22.1 pm)/En/i (3.35 pm)] 
and [WISE Fw 3 (11.6 pm)/En/i (3.35 pm)] as functions of 
[2MASS Fj (1.25 pm)/Ejf^ (2.17 pm)] for 79 Algols having po¬ 
sitional matches with both 2MASS and WISE within less than 
3" (31 objects have matches better than 0712 and 15 better than 
0706). Spurious detections due to a diffraction spike from a close 
bright star, the persistence from a bright star scanned immedi¬ 
ately before the target, halo scattering from a close source or 
optical ghost have all been discarded. We retained systems hav¬ 
ing a signal-to-noise ratio (SNR) in band W4 higher than 3.2. 
In addition, we note that the variability labels range from ‘not 
variable’ (about 45 objects; this number depends on the band 
considered) to ‘high probability of being true variable’ (34 ob¬ 
jects). This variability observed by WISE could be attributed to 

http://cdsarc.u-strasbg.fr/viz-bin/Cat?II/125 
bttp://cdsarc.u-strasbg.fr/viz-bin/Cat?II/156A 
''' bttp: //cdsarc. u- strasbg. fr/viz-bin/Cat?II/246 
bttp://cdsarc.u-strasbg.fr/viz-bin/Cat?II/311 
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Fig. 13. Top panel: WISE Fy/ilFwi against 2MASS J/Kj flux ratios for Algols from the catalogue of Budding et al. (2004), and for W Ser 
from the list of Gudel & Elias (1996). The symbiotic Algol SS Lep and a few Be and B[e] systems have been considered as well (as discussed in 
Sect. 4.1). The black filled triangles (linked by the dashed line) are black bodies of given temperatures (labels in Kelvin). The red filled squares 
represent the simulations (STAR: input spectrum). For clarity, only extreme models are labelled, we refer to Table 5 for unlabelled models. The 
filled circles correspond to Algols with a SNR > 3.2. The size of the circle is proportional to the SNR for band W4 (except for SS Lep which has 
a SNR of 962 and for which we limit the size of the symbol) and the symbol colour relates to the for the same band. Bottom panel: Same as 
top panel but for FwilF^i with SNR and^^ relative to IV3. 


an outflow, but it could also simply result from the eclipse of 
the donor by the gainer (in our model, the donor is three times 
brighter than the gainer in the W4 band). A more precise anal¬ 
ysis of variability with WISE data is beyond the scope of this 
study. 

Figure 13 includes observed systems (filled circles) and our 
models (red squares). The circle colour codes the (which 
measures the deviation from the point-spread function-PSF; it 
is larger than 1 for extended sources) for bands W4 and W3. The 
circle size is proportional to the SNR in bands Wh and WA-. The 
J and Ks bands are not affected by dust emission (see Fig. 9), 
and only depend on the stellar temperatures and on the covering 
factor. Since our simulations are based on one specific Algol 
system at three different stages of its evolution (see details in 
Table 1), the range in FjfF^^ values covered by the models is 
limited compared to observations. In the adopted Algol model 
system, the donor star dominates the spectrum at longer wave¬ 
lengths, beyond the J band (Figs. 9 and 12). 

On the other hand, FwaIF^i directly traces the presence 
of IR emission from dust grains (Fig. 9), and Fw?,/F^i probes 
the silicate emission around 12 pm. The models distribute 
along three horizontal sequences in the upper panel of Fig. 13, 
dust-free models corresponding to the lowest F^aIFwi value 
(close to the stellar locus), followed by low-dust models at an 
intermediate FwaIFwx value, and finally the high-dust model 
’C_hd_0.71_inter’ falling above the limit of the figure. 

Fig. 13 reveals that several Algols and a few W Ser and (3 Lyr 
in our sample have FwaIFwx ratios falling above the black-body 
sequence, thus hinting at the presence of dust emission. The 
same holds true for the FwaIFwi ratio, which has not been dis¬ 
played though because the 1T2 flux may reach (close to) satura¬ 
tion for bright systems. For four Algols (CZ Vel, UZ Nor, DM 
Pup and V509 Mon), one p Lyr (RY Set), and the symbiotic Al¬ 


gol SS Lep, FwaIFwi is close to or larger than 1, well above 
our model predictions for a dust-to-gas ratio of O = 2 x 10“^. 
The SNRs of the WISE measurements for these systems are 
among the highest in our sample (see Table 6). The WISE and 
2MASS data for these systems'® and others with moderate ex¬ 
cesses are collected in Table 6. Among these, not much is known 
(apart from their eclipsing light curve) for DM Pup, V509 Mon, 
CZ Vel, UZ Nor, V442 Cas, and DH Her, but given the interest¬ 
ing nature of all the other stars picked up by the F^nlF^x crite¬ 
rion, a dedicated study of this small list of not well-studied stars 
will certainly pay off. RY Set is a very interesting contact binary 
consisting of two seemingly O-type supergiants (but the 11 d or¬ 
bital period makes it unlikely that the components are really as 
extended as supergiants are), with dust forming in an ionising 
radiation field, and with a radiation-driven nebular wind with a 
kinematic age of 120 yr (e.g., Smith et al. 2001). 

Finally, three more systems have Fw 4 /Fwi ratios close to 
the value expected for O = 2 x 10 the Be systems V696 Mon 
(= HR 2142), 4> Per and SX Aur. SX Aur shows an extended 
shell in the W4 band as can be seen in Fig. 14. This is an in¬ 
teresting result in view of the controversy regarding the exact 
evolutionary status of that system: from the light curve analysis, 
Chambliss &. Leung (1979) obtained a semi-detached configu¬ 
ration with the secondary filling its Roche lobe, and endorse an 
earlier suggestion by Stothers (1973) that the system is under¬ 
going case A mass transfer and that mass reversal has occurred. 
On the contrary, based on a new light curve analysis, Linnell et 

We do not confirm thus the conclusion of Taranova & Shenavrin 
(1997) who seemingly found evidence for dust around the W Ser sys¬ 
tem RX Cas, though at shorter wavelengths (in the LM bands at 3 - 5 
pm). For this system, we found no significant infrared excesses from 
WISE (T’m' 4 /Th/i = 0.056). Therefore, the excess found by Taranova & 
Shenavrin (1997) in the LM bands cannot come from dust. 
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Table 6. 2MASS and WISE magnitudes for the 13 systems deviating the most from the hlack-hody law (see Fig. 13). 


Name 

Poi-b 

(Sp.l) 

J 

H 

K, 

W1 

W2 

W3 

W4 


(d) 


(mag.) (mag.) 

(mag.) 

(mag.) 

(mag.) 

(mag.) 

(mag.) 


q 

[Sp.2] 

vj (mas) 



[SNRwi ]Cfwi) 

[SNRw2](Tw2) 

[SNRw3](Tw3) 

[SNRw4](Tw4) 


5.1927d) 

(B3)(3) 




Algols 




CZ Vel 

9.441 

9.192 

9.032 

8.562 

8.485 

7.814 

3.669 


0.50® 

[B6] 

Unk. 



[50](1.75) 

[52.6](1.79) 

[48.2](5.46) 

[42.9](35.4) 

UZ Nor 

3.1960® 

(A1)(3) 

10.262 10.079 

9.982 

9.873 

9.817 

8.202 

5.937 


0.18® 

[KOIV] 

Unk. 



[48.3](3.27) 

[57.4](2.82) 

[40.8](5.7) 

[22.9](1.33) 

V509 Mon 

4.9174® 

(B7)(3) 

12.595 12.171 

11.987 

11.868 

11.905 

11.515 

8.38 


0.32® 

[G4IV] 

Unk. 



[46.3](1.23) 

[43.9](1.37) 

[6.5](1.14) 

[4.7](0.992) 

DM Pup 

3.567''3) 

(B5)(3) 

12.761 

12.569 

12.488 

12.057 

12.035 

11.215 

8.868 


Unk. 

[A2.5] 

Unk. 



[46.1](1.38) 

[42.6](1.61) 

[9.6](1.00) 

[3.2](1.09) 

V442 Cas 

3.5922(3) 

(A7)(3) 

11.729 11.516 

11.407 

11.286 

11.271 

10.961 

8.76 


0.18® 

[K2IV] 

Unk. 



[49.6](1.36) 

[52.2](1.42) 

[14.9](0.995) 

[4.2](0.994) 

DHHer 

Unk. 

(A5)(6) 

10.054 

9.997 

9.956 

10.134 

10.168 

10.189 

8.449 


Unk. 

[Unk.] 

Unk. 



[47.9](5.90) 

[53.6](6.06) 

[26.9](1.12) 

[5.4](0.979) 


11.12® 

(09.7Ibpe)(i')) 




ySLyr 




RYSct 

6.245 

5.854 

5.463 

4.955 

3.895 

0.014 

-0.922 


0.23("( 

[06.5 I] 

1.59 



[14.2](12.1) 

[19.4](51.0) 

[21.3](0.550) 

[83.96](2.21) 


18.6('4) 

(A7I)(3) 




WSer 




V367 Cyg 

5.202 

5.006 

4.828 

4.72 

4.038 

3.781 

2.676 


0.825'® 

[Unk.] 

1.74 



[13.2](0.434) 

[20.9](2.34) 

[70.9](1.27) 

[31.4](1.24) 

WSer 

14.154® 

(B emb.)(i3) 

7.416 

6.109 

5.558 

5.106 

4.560 

4.038 

3.379 


0.64'3) 

[F5III] 

0.70 



[13.6](0.737) 

[18.2](3.65) 

[69](1.72) 

[31.3](1.36) 


260.3(3) 

(A1V)(3) 



Symbiotic Algol 




SS Lep 

2.094 

2.105 

1.671 

3.455 

1.536 

-1.307 

-2.117 


0.47(3) 

[M6II] 

3.59 



[39](186) 

[64.3](195) 

[3](0.0094) 

[962] (45.2) 

HR 2142 

~80('3) 

(B2IVne)(i3) 

5.088 

5.068 

4.781 

Be 

4.828 

3.756 

3.31 

2.592 


Unk. 

[Unk.] 

2.48 



[12.1](15.4) 

[16](25.5) 

[91.1](7.12) 

[58.9](1.96) 

TT Aqr 

84.1(16) 

(BlVe)('i) 

5.305 

5.365 

5.351 

4.573 

3.767 

4.056 

3.58 


0.16(16) 

[A-F](16) 

4.17 



[13.8](1.22) 

[16](2.10) 

[73.8](1.68) 

[41.7](1.02) 

CX Dra 

6.696(13) 

(B2.5Ve)(i3) 

5.870 

5.805 

5.486 

5.457 

5.140 

5.128 

4.854 


0.23(16) 

[F5III](13) 

2.52 



[17.4](1.27) 

[32.9](2.07) 

[77.9](4.48) 

[49](1.11) 

(p Per 

127(13) 

(B2e)(i3) 

4.049 

3.955 

3.709 

3.68 

2.774 

2.639 

1.869 


0.158(3) 

[sdO] 

4.54 



[8.7](0.247) 

[12.2](1.59) 

[53.5](2.03) 

[65.4](0.737) 

SX Aur 

1.2l('i) 

(B2)('i) 

8.431 

8.471 

8.541 

8.594 

8.655 

8.636 

7.047 


0.61(3) 

[B5] 

1.39 



[47.9](2.97) 

[55.6](3.53) 

[42](1.29) 

[13.6](1.34) 


Notes. ‘Unk.’ refers to unknown values, ‘emb.’ to an embedded star where the spectral type is only roughly guessed, q is the mass ratio defined as 

^ donor/gainer- 

References. Zasche 2009; Zasche 2011; Budding et al. 2004 and references therein; Linnell et al. 1988; Blind et al. 2011; *** Kreiner 
2004; <’> Isles 1985; ™ Linnell et al. 1988; Kreiner 2004; Walborn 1982; '*'> Grundstrom et al. 2007; <'2) Plavec 1989; Hegedus 1988; 
Zola & Ogloza 2001; Hummel & Stefl 2001; Bjorkman et al. 2002; Waters et al. 1991; Berdyugin & Piirola 2002. 


al. (1988) concluded that the system is made of two unevolved B 
stars, the primary being however on the verge of filling its Roche 
lobe. Our discovery of an extended infrared shell around SX Aur 
(Fig. 14) seems to suggest that some mass was lost from the sys¬ 
tem in a recent past, even though mass transfer is not currently 
occurring. This result confirms the interest of considering Be 
stars as well when looking for evidence of systemic mass loss 
from Algols (see Sect. 4.1). Concerning HR 2142, it must be 


noted that the target lies close to a star forming region (Mono- 
ceros R2 at 56') and it is likely that the infrared emission does 
not originate from the binary system itself but rather from the 
rich sky background. 

Among the systems with the highest FwaIFwi ratios, we 
note that CZ Vel is the only genuine Algol system with an ex¬ 
tended nebula detected by WISE (Fig. 14), as could be inferred 
by its highx^ 4 - Its largest angular size is about 180". Unfor- 
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Fig. 14. WISE band-4 images for the genuine Algol CZ Vel (left) and the Be star SX Aur (right). The images are centred on the objects, with 
North up and East left. The flux scale goes from the background level to background plus 150 times the background scatter cr, and the first thin 
(green) contour line corresponds to the background level plus 3 times its fluctuation, the next intermediate (blue) contour curve to 5 cr and the 
inermost thick black contour line to 10 cr. 


tunately, the parallax of CZ Vel is unknown, as is therefore the 
linear size of the nebula. The nebula is very asymmetric and 
does not resemble our predictions (Fig. 12). 

There is something else to remark from Fig. 13. Systems 
like W Ser or V930 Oph (having the lowest FjIFk, - 0.4 and 
0.71, respectively) exhibit very low equivalent black-body tem¬ 
peratures {T < 3000 K), in fact well below the components’ ef¬ 
fective temperatures. These low temperatures thus indicate that 
the system light is, at those infrared wavelengths, dominated by 
a very cool object, possibly an accretion disc which obscures 
the gainer star, as in W Ser. The very short period {P = 1^4; 
Maikov 1993) of V930 Oph does not support, however, the pres¬ 
ence of an accretion disc; direct impact of the accreting stream 
is expected instead (Lubow & Shu 1975). The low-temperature 
black body detected is therefore likely to originate from a sub¬ 
stantial amount of cool dust far from the system, indicative of 
systemic mass loss. 

Figure 13 shows that a typical ISM dust-to-gas ratio (O = 
4 X 10“^, model C_hd_0.71) produces F^aIFwi ratios larger 
than those observed with WISE in Algols with an infrared ex¬ 
cess. Thus, observations favour a smaller dust-to-gas ratio in 
the stellar nebula than in the ISM. This conclusion is consis¬ 
tent with dust being destroyed by high-energy photons from 
both the B star and the hotspot. It is likely for the same reason 
that so few genuine Algols exhibit infrared excesses caused by 
dust. Moreover, in the framework of the hotspot mechanism (Pa¬ 
per I), systemic mass loss occurs in the early stage of mass trans¬ 
fer. Since the circumstellar shell dissipates very quickly (within 
300 years) after the cessation of the non-conservative evolution, 
this is another reason why infrared excesses are expected to be 
rare around genuine Algols. 

Similarly, W Ser systems rarely show the signature of dust 
emission. This is compatible with our finding that the presence 
of highly ionised species (as seen in W Ser systems) is exclusive 
of the presence of dust (see Sect. 3.3). 

4.3. Constraint on Stromgren photometry and weakly ionised 
emission lines 

As mentioned in Sect. 3.2.1, the presence of a strong Balmer 
continuum in emission and the forest of metallic-line emission 
in the emergent spectrum are not typical of Algols. Stromgren 


photometry, through its c\ index (Sect. 3.2.5 and Fig. 10), offers 
an easy way to check for the presence of Balmer emission con¬ 
tinuum. Usually, as concluded before, IR excesses and emission 
lines or continuum are mutually exclusive. There are three inter¬ 
esting exceptions, though: the (3 Lyr system RY Set, consisting 
of two stars with O-type spectra (Smith et al. 2001) and dust 
formation in the circumstellar nebula despite ionising radiation 
(Gehrz et al. 1995). We confirm the latter authors’ finding, since 
RY Set is the system with the largest FwaIF^i ratio (Table 6) 
and at the same time the smallest ci_o index (Table 7), indicative 
of a strong Balmer continuum in emission. The system clearly 
has all features indicative of systematic mass loss. The Be stars 
(f) Per and n Aqr also present a negative ci_o, indicative of a strong 
Balmer continuum. This property thus reveals once more the in¬ 
terest of Be binary systems in the search for Algol-like systems 
with systemic mass loss. 

In Table 7, there are 3 other objects with negative ci,o values, 
namely the W Ser stars AX Mon and V644 Mon, and the Be star 
HR 2142. All the other systems are located in a region of the m\ Q 
- ci_o diagram not typical of Balmer-continuum emission. On 
the contrary, the W Ser star V367 Cyg has a large ci_o and devi¬ 
ates considerably from the main sequence locus. Young & Sny¬ 
der (1982) presented a dedicated study off uvby photometry of 
V367 Cyg (along with a few other W Ser systems) and found the 
same peculiar Stromgren colours as those reported here, which 
imply either excess flux longward of 400 nm or a strong absorp¬ 
tion in the u band, due for instance to improperly corrected in¬ 
terstellar absorption. The latter interpretation is made likely by 
the fact that V367 Cyg is a remote object in the Cygnus dust 
clouds. Therefore, Young (& Snyder (1982) conclude that “there 
is reason to question the reddening correction” in such a peculiar 
environment. Of course, one cannot firmly exclude an alterna¬ 
tive reason related to the physics of the binary system, despite 
the fact that the (mi o, ci_o) colours of V367 Cyg are not typical 
of the other W Ser systems studied here and by Young & Snyder 
(1982). 

5. Conclusions 

We performed photoionisation and radiative-transfer simulations 
of the material surrounding an Algol system that is loosing mass 
via the hotspot mechanism. We calculated models with various 
configurations (different hotspot temperatures, with or without 
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Table 7. Stromgren mj o and ci,o for our sample of W Ser and Be stars. 


Name 

Type 

mi,o (mag.) 

ci,o (mag.) 

Ref. 

RYSet 

/3Lyr 

-0.0218 

-0.3002 

(4) 

phyr 

/3Lyr 

0.0697 

0.32 

(12) 



0.07006 

0.3186 

(2) 



0.06518 

0.3278 

(6) 

AUMon 

WSer 

0.05866 

0.3186 

(1) 



0.04928 

0.2748 

(5) 



0.0599 

0.144 

(10) 



0.06386 

0.3976 

(7) 

KXAnd 

WSer 

0.01472 

0.1662 

(7) 

RZ Set 

WSer 

-0.00166 

0.1054 

(1) 



-0.01142 

0.0768 

(9) 



0.01804 

0.1594 

(6) 

UX Mon 

WSer 

0.14302 

0.3742 

(10) 

V356 Sgr 

WSer 

0.0435 

0.527 

(9) 

V367 Cyg 

WSer 

0.03818 

1.4508 

(1) 



0.03744 

1.4494 

(7) 



0.0383 

1.445 

(8) 



0.03424 

1.4944 

(6) 

WCru 

WSer 

0.41384 

0.4614 

(1) 



0.41672 

0.4352 

(14) 



0.4062 

0.538 

(9) 

WSer 

WSer 

0.14036 

0.3076 

(9) 

V395 Aur 

Be 

0.15024 

0.2954 

(6) 

V617 Aur 

Be 

0.14018 

0.2518 

(7) 

V644 Mon 

Be 

-0.01952 

-0.0202 

(10) 

V1914 Cyg 

Be 

0.17682 

0.1162 

(13) 

UCrB 

Be 

0.07786 

0.4606 

(6) 

AX Mon 

Be 

-0.02574 

-0.1174 

(6) 

(p Per 

Be 

0.0336 

-0.095 

(4) 

CX Dra 

Be 

0.06112 

0.2262 

(15) 

;r Aqr 

Be 

0.00188 

-0.1892 

(1) 



-0.00296 

-0.2096 

(2) 



0.0038 

-0.179 

(3) 

HR 2142 

Be 

0.02558 

-0.0042 

(1) 



0.0236 

0.017 

(2) 



0.02648 

-0.0152 

(3) 


References. Hauck & Mermilliod 1998; Grpnhech & Olsen 1976; 

Crawford et al. 1971; Westin 1982; Knude 1992; Cameron 

1966; Olsen 1983; Perry 1969; Wolf & Kern 1983; Kilkenny 

et al. 1985; Reglero et al. 1987; Hauck & Mermilliod 1980; 
Olsen 1993; Gray & Olsen 1991; Crawford et al. 1971. 

dust, different opening angles for the disc-like outflow), and find 
that the resulting IR excesses produced by dust in the outer parts 
of the outflow are consistent with WISE observations for a small 
fraction of Algols, especially the systems CZ Vel, UZ Nor and 
V509 Mon. We also report infrared excess in the Be system 
SX Aur, in the p Lyrae system RY Set and in the symbiotic Algol 
SS Lep (already known for the latter two). As expected, most 
genuine Algols do not exhibit strong infrared excesses since the 
systemic mass loss occurred at an earlier phase and the mate¬ 
rial has since dissipated. Our models also predict strong Balmer 
continua that are not observed in Algols but are visible in some 
Be systems. 

The outflowing gas can significantly reduce the observed UV 
continuum-flux level. We suspect this absorption to quickly dis¬ 
appear after the end of the non-conservative phase due to the 


high outflow speed of the line-driven wind and to show some 
phase dependency as the matter is launched from the hotspot. 
This variability may help to identify systems undergoing sys¬ 
temic mass loss due to the hotspot mechanism. 

Our simulations indicate that the geometry of the outflow has 
an impact on the emergent SEDs, especially in the IR range for 
dusty models, where the fluxes decrease with increasing opening 
angle of the outflow. 

In the SEDs computed from our models with a hotspot tem¬ 
perature Ths ^ 35 000 K, we find strong emission lines of highly 
ionised species such as Si iv which arise close to the stellar sur¬ 
face, inside the Roche lobe of the gainer star, and are typical of 
W Ser systems. We emphasise the need for further observations 
of the extended family of Algols, which includes DPVs, sym¬ 
biotic Algols, Be and B[e] binaries, W Ser and p Lyr systems. 
These types of objects can be related to different steps in the 
evolution of Algols or to different initial configurations (masses, 
periods). The search for signatures of a non-conservative evo¬ 
lution should not focus specifically on Algols, since they are 
older Class III systems (Paper I), but more generally on Algol- 
related systems where systemic mass loss can potentially occur. 
We also hope that this first study of the hotspot-outflow for¬ 
mation will motivate a more precise analysis of the stream/star 
or stream/disc interaction, preferentially by use of 3D mod¬ 
els including radiation-hydrodynamics (see Bisikalo et al. 2000; 
Nazarenko & Glazunova 2013, for p Lyr). 
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Appendix A: Estimation of the opening angie of the 
outflow 

The geometry of the escaping stream may be complex, forming 
a disc-like geometry, and possibly creating a spiral-like structure 
close to the star (see Sect. 2.2). The geometrical thickness of the 
disc structure is determined by the opening angle of the outflow. 
Since the outflow is driven by the radiative flux of momentum 
Pv, we can estimate the geometrical thickness of the outflow by 
use of the dependence of Pv on the direction given by (an¬ 
gle measured with respect to the direction perpendicular to the 
hotspot surface; see Eig. A.l). 

The total momentum flux, normal to an element of surface 
that is parallel to the hotspot, is given by (Eig. A. 1) 

Pv(6n = 0) = — flh r cos^(0)dQ, (A.l) 

^ .7hemisphere 

where ly is the specific intensity (assumed to be isotropic) orig¬ 
inating from the hotspot, fih is the unit vector normal to the 
hotspot, c is the speed of light and dQ = sin ff d0 d(p. 

To calculate the momentum flux toward an arbitrary angle 
0n < 7r/2, we introduce the reference frame (x,'y,'z), where the Y 
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Fig. A.l. Geometry of the hotspot emission, iih is the normal to the 
surface of the hotspot and iijjyi is the normal of the surface d^, which 
forms an angle with fli,. 6flow/2 (see text) is equal to when the mo¬ 
mentum flux integrated in azimuthal angle between 0 and equals 68% 
of the momentum flux integrated over the total hemisphere (between 0 
and |). 


axis lies in the plane of the hotspot, the ^ axis lies in the plane 
d^, andTis collinear with fldja, the normal to d^. Note that 0 is 
measured with respect to The general expression of Eq. A. 1 
for arbitrary angles 0„ reads 

Pv(0n) = —fldja r cos 6 |cos 9\ 0(n ■ flh) dQ, (A.2) 

../hemisphere 

where n is the unit vector pointing in the direction {9, (p). The 
absolute value in the second factor of the integrand ensures that 
radiation traversing the surface d^ from above adds a negative 
contribution to the total momentum flux. The Heavyside step 
function 0(n • flh) accounts for the fact that the hotspot is only 
radiating upwards. 

For a fixed angle 0n, the vector flh normal to the hotspot has 
coordinates (0, - sin6>n, cos so that an evaluation of the dot 
product gives the limits of integration in (p: 


arcsin [ 1 / (tan 9 tan 0n)] 

9 > 7r/2 - 9a 

, (A.3) 

0 

otherwise 

7T - arcsin [l/(tan6/tan0n)] 

9 > njl - 9a 

, (A.4) 

2n 

otherwise 


where the latter expression for (pum .2 results from the symmetry 
of the problem (the angle ip is 0 above the T axis). With these 
limits, Eq. A.l can be expressed as 


Pv(0n) 



(A.5) 

(A.6) 


(A.7) 


where the first integral (Eq. A.5) corresponds to the hemisphere 
over the surface element d^ in Fig. A.l and the second integral 



Fig. B.l. Representation of the effective surface of the hotspot over a 
full orbital period. See text for description. 


(Eq. A.6) corresponds to the two green parts. The first of the 
latter is subtracted because it lies below the hotspot surface and 
does not contribute, the second one has a negative contribution 
to the momentum flux because it is below the surface element 
d^. These two regions are identical by symmetry and allow us 
to use the same limits of integration in (p and 9. 

Our estimate of the opening angle 6/flow of the outflow is 
based on the integral of the momentum flux Pv over 0n- We de¬ 
fine it as the angle that corresponds to a fraction of 0.68 (arbi¬ 
trarily set; standard 1 cr as for a Gaussian distribution) of the total 
integral, i. e., 

^flaow/2 prrl2 

P,(e„)d0„ =0.68 Pv(0n)de„. (A.8) 

Jo Jo 

This gives an opening angle of 6/flow ~ 90°. The ‘effective’ solid 
angle in which the disc-like outflow spreads is given by 

p2K ^7r/2+Sjow/2 

A-nfc - I dip I d0sin0 = 47rsin(6/flow/2) ~ 2.847rsr, 

Jo J;r/2-efl„w/2 

(A.9) 

where /c is the covering factor that is defined as the fraction of 
the spherical domain that is covered by gas (Sect. 2.1). This is 
different from the geometry adopted in Paper I, where the out¬ 
flow was arbitrarily assumed to be confined in a quarter of a 
sphere (solid angle equal to n). 


Appendix B: Calculation of the effective flux 

In our simulations, the effective radiative flux Eetj is given by 

T’eff = Eg + (B.l) 

-^eft.hs 

where Eg = and Ehs = Ehs/.^hs are the gainer star and 

hotspot fluxes, respectively, JI* is the surface area of the star, 
and Jlfls is the surface area of the hotspot, i. e. the surface of a 
calotte with the radius of the base equivalent to the radius of the 
impacting stream, 

JIhs = 27rE, , (B.2) 
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Wavelength (Angstrom) 



Fig. C.l. Top: Continuum in the UV hand for W Ser observed with 
lUE. Black: SWP47439 observation made on 1993-04-07 (20:02:16) 
at phase O = 0.84; red: SWP47457 observation made on 1993-04-10 
(21:40:17) at phase 4> = 0.06. Bottom: Difference of the two spectra. 


Fig. C.2. HST/GHRS spectra of the [Siiv] 1393.76 and 1402.77 A 
lines of W Ser at two distinct phases (black: z01u5108t observation 
made on 1991-07-12 (05:32:25) at phase O = 0.00, primary eclipse; 
red: z01u0308t observation made on 1991-07-19 (03:29:13) at phase 
O = 0.5, secondary eclipse). 


where R, and rhs are the star and hotspot radius, ^efi.hs is the area 
of the strip at the equator of the star of height 2rhs (Fig. B. 1), 

^eff,hs = 47r/?J - 2 X 2nRt{Rt - rj,.,), (B.3) 


over which the hotspot flux is averaged. Therefore, the effective 
flux reads 


FeS - Fg H- 


2rhs 




(B.4) 


The net flux from the hotspot is thus reduced by the factor 


/?. 

2rhs 




(B.5) 


In our simulations, the source of radiation emits isotropically, so 
that the total effective luminosity will be 

Teff — Tg -I- ~ Tg -H (B.6) 

rhs 

where Fgff hs ~ (?^hs/?^eff,hs)Fhs. 


at phase 0.1 during the primary eclipse. We note that at phase 
0.1, the absorption in the longer wavelength region of the spectra 
(A > 1400 A) and the Lyman a emission are stronger than during 
the secondary eclipse, possibly indicative of a hot ionised region. 

Figure C.2 presents the profiles of the [Siiv] 1393.76 and 
1402.77 A lines as observed with HST/GHRS for W Ser at two 
different phases: primary eclipse (d) = 0.0), and secondary 
eclipse (<t) = 0.5). The same ephemeris has been used. We 
clearly see an excess of flux in the blue wing at phase <t) = 0.5. 
The total extent of the wings of both lines reaches 600 km s ', 
which is comparable with the escape velocity at the surface 
of the star (643 kms“* considering a main sequence gainer 
star of 1.51 Mq according to Budding et al. 2004, and follow¬ 
ing a mass-radius relation for main-sequence stars of the form 
R oc m"-*) and with the outflow terminal velocity in our model 
(850 km s '). In their attempt to model this line variability, Wei- 
land et al. (1995) found two components, a broad one likely due 
to the boundary layer of an accretion disc and a narrow one pos¬ 
sibly due to the formation of a hotspot. Further investigation, 
especially at different phases, may shed new light on the origin 
of the [Si iv] emission. 


Appendix C: UV continuum and iine variabiiity in 
WSer 

W Ser systems are known to show a strong line variability in 
the UV (especially for highly ionised species; see the review by 
Plavec 1982). In this appendix, we investigate, for the proto¬ 
typical system W Ser, the observed UV spectra in the search for 
variations in the continuum as presented in Sect. 3.2. We also 
demonstrate the variability of [Si iv] lines for this object. 

Figure C.l shows low-dispersion, large aperture spectra of 
W Ser in the UV as observed by the International Ultraviolet 
Explorer (lUE) at two distinct phases (these phases are estimated 
thanks to the improved ephemeris of Piirola et al. 2005). The 
two spectra have been taken only 3 days apart limiting the effect 
of intrinsic variability. The red spectrum is taken between the 
secondary and the primary eclipse at phase 0.8 and the black one 


Appendix D: Chemical composition of the outflow 

Table D.l presents the chemical composition of the donor star 
atmosphere that will be transferred by RLOE onto the accretor’s 
surface and eventually ejected by the hotspot. 
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Table D.l. Chemical composition of the outflowing material for the three models computed. 


Element (i) 

(Model A) 

(Model B ) 

(Model C) 

(Solar ref.) 

H 

6.9941 X 10 ‘ 

6.9902 X 10' 

6.9810 X 10' 

7.3732 X 

10 -' 

He 

2.8059 X 10-' 

2.8074 X 10-' 

2.8134 X 10-' 

2.4917 X 

10 -' 

Li 

6.7290 X lO-'*' 

3.1061 X 10-'3 

1.0826 X 10-'3 

1.0358 X 

10-3 

C 

3.4675 X 10-3 

2.2462 X 10-3 

1.0944 X 10-4 

2.1526 X 

10-3 

N 

1.0632 X 10-3 

2.5270 X 10-3 

4.9937 X 10-3 

8.7194 X 

10-4 

O 

9.6519 X 10-3 

9.6486 X 10-3 

9.6333 X 10-3 

5.7348 X 

10-3 

F 

5.6122 X 10-3 

5.6411 X 10-3 

5.6984 X 10-3 

4.1970 X 

10-3 

Ne 

1.9697 X 10-3 

1.9696 X 10-3 

1.9575 X 10-3 

1.4762 X 

10-3 

Na 

4.0004 X 10-3 

4.0130 X 10-3 

5.2729 X 10-3 

3.5989 X 

10-3 

Mg 

7.5221 X 10-'' 

7.5221 X 10-4 

7.5220 X 10-4 

6.1694 X 

10-4 

Al 

6.4825 X 10-3 

6.4825 X 10-3 

6.4825 X 10-3 

5.8225 X 

10-3 

Si 

8.1136 X 10-4 

8.1136 X 10-4 

8.1136 X 10-4 

7.1290 X 

10-4 

P 

7.1102 X 10-3 

7.1102 X 10-3 

7.1102 X 10-3 

7.2504 X 

10-3 

S 

4.3015 X 10-4 

4.3015 X 10-4 

4.3015 X 10-4 

4.3159 X 

10-4 

Cl 

heavy 

2.3983 X 10-3 
1.7395 X 10-3 

2.3983 X 10-3 
1.7395 X 10-3 

2.3983 X 10-3 
1.7395 X 10-3 

4.9535 X 

10-3 


Notes. Contributions from different isotopes (when available) are added up. For the Cloudy simulations we use the thirty lightest elements. For 
those not appearing in the table we adopt solar abundances. 

References. Solar abundances: Grevesse & Sauval (1998), except N, Ne, Mg, Si, and Fe (Flolweger 2001), O (Allende Prieto et al. 2001), and C 
(Allende Prieto et al. 2002). 
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